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1 Objectives 

The main objectives of the analytical part of this investigation are to study the fluid 
flow phenomena together with the thermal effects on drops levitated in an acoustic field. 
To a large extent, experimentation on ground requires a strong acoustic field that has a 
significant interference with other thermal-fluid effects. While most of the work has been 
directed towards particles in strong acoustic fields to overcome gravity, some results for 
microgravity have been obtained. 

One of the objectives was to obtain the thermocapillary flow in a spot-heated drop, and 
set up a model for the prediction of thermophysical properties. In addition, for acoustically 
levitated particles, a clear understanding of the underlying fluid mechanics was required. 
Also, the interaction of acoustics with steady and pulsating thermal stimuli was required 
to be analyzed. 

The experimental part of the work was funded through JPL, and has been reported 
separately. 


2 Accomplishments 

The analytical developments have led to several new and interesting results pertaining to 
heat transfer and fluid mechanics. A summary of each set of results is given here. 

2.1 Thermocapillary Flow in a Spot-Heated Drop 

In microgravity situations, the low Marangoni number flow for a spot-heated liquid drop 
has been analyzed. As a first approximation, heat flow treated as pure conduction. The 
liquid region is considered to be fully transient while the surrounding gas medium is 
regarded as quasisteady. From the analytical solution for the temperature distribution, the 
thermocapillary force on the interface was calculated and applied as an interface condition 
for Stokes flow in both the media. This led to the complete flow streamlines in the drop, 
as well as an expression for the migration velocity. The results have been presented for a 
liqiod-to-solid conductivity ratio of 100 and Fourier numbers ranging from 0.1 to 10. The 
inference of the the thermal diffusivity form the calculations have been discussed in the 
paper which has been accepted for publication in Microgravity Science and Technology , 
Vol. X (1997). 

2.2 Acoustic Streaming Around a Sphere 

For a particle levitated in an acoustic field, its position is determined, to some extent, by 
the compressibility properties of the particle phase relative to the surrounding medium. 
For example, liquid particles in a gas tend to position themselves at the velocity antinode 
while gas bubbles equilibrate at the velocity node. The effect of a gravity-type body 
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force may offset the equilibrium position to a point somewhere between the node and the 
antinode. 

The fluid mechanics for the case of a particle at the velocity antinode can be derived 
directly from Riley’s (1966) work for a vibrating sphere in an otherwise quiescent fluid. 
However, the flow about the velocity node has not been available and a detailed analysis 
for that situation has been carried out. In the present development, the perturbation 
procedure of Riley (1966) is employed to derive the flow field for the situation when a 
spherical particle is positioned at the velocity node. As in Riley’s solution which applies 
to a sphere at the velocity antinode, it is found that there is a thin shear-wave region 
adjacent to the spherical boundary. However, for this thin Stokes layer, the streamlines 
are not closed but join with the outer flow. Therefore, the shear- wave layer does not cover 
the entire sphere as in the previous case, but lies mostly around the equatorial region of 
the sphere. The equatorial belt lies in the region 55° < 0 < 125°. 

Besides providing detailed knowledge about the flow field around the particle at the 
velocity node, the present analysis is very useful for developing the solution when the 
particle between the node and the antinode. This will be a combination of the two 
solutions along with some nonlinear interaction terms. 

The velocity-node solution is to be submitted for publication. The work has been 
presented as a part of Ms. Hong Zhao’s Ph.D. dissertation proposal. 

2.3 Internal Circulation in a Drop in an Acoustic Field 

The internal flow in a drop at the antinode of a standing wave has been investigated. 
This was done through the application of the stress continuity condition at the liquid-gas 
interface. To the leading order of calculation, the internal flow field was found to be very 
weak. At the next higher order, steady internal flows are predicted. There is, however, 
an important effect on the recirculating Stokes layer which vanishes when 

|M| < (§) >/2[2 + 5 (£/,*)]> 

where fi is the liquid viscosity, p is the exterior gas-phase viscosity, and M is the dimen- 
sionless frequency parameter for the gas phase, defined by 

, iua 2 p 

M = . 

ft 

This is a very interesting new development which needs detailed confirmation with exper- 
iments. While it agrees with many expriments with levitated drops where no recirculating 
layer has been found, a new set of experiments for specifically testing the theory need to 
be cerried out. 

The work has also been presented as a part of Ms. Hong Zhao’s Ph.D. dissertation 
proposal. 
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2.4 Thermoacoustic Streaming from a Sphere 

This work consists of a study of the thermal effects of the steady streaming motion induced 
by a solid sphere. The thermal field is generated by an oscillating boundary condition on 
the surface of the sphere. The interesting feature here is that when the thermal frequency 
is ‘in tune’ with the acoustic field, the convective heat transfer term remains nonzero when 
the energy equation is averaged over time. This leads to a net flow of energy through the 
equatorial plane of the sphere. 

A singular perturbation analysis has been carried out for small values of the streaming 
Reynolds number (U^/Ru 1), along with \M\ = d 2 uj/u » 1. There is equal heating 
and cooling of the regions above and below this plane, respectively. A result for the net 
flow across the equatorial plane of the sphere has been obtained. 

The above case of thermoacoustic streaming has been extended to the case of pulsating 
spot heating. Again, the frequency of oscillation is in tune with the acoustic field and 
the time-averaged energy equation describes nonzero convective transport. The thermal 
energy deposited is divided between the solid sphere and the surrounding gas medium 
on the basis of the respective thermal transport characteristics. The calculations have 
provided an analytical expression for the transient thermal field. The analysis tor long 
times (large Fourier number) predicts steady buildup of the time-averaged temperature 
in the sphere. While experimentation with thermoacoustics has shown somewhat similar 
behavior, there is some skepticism about the present prediction which could perhaps be 
an artifact of approximation. 

A paper (Gopinath & Sadhal, 1994) resulting from this work has been published in 
the Proceedings of the 1994 International Heat and Mass Transfer Conference (Brighton, 
UK). 

3 Microgravity Sessions in ASME Conferences 

The Principal Investigator has taken a role in organizing technical sessions the National 
Heat Transfer Conference and the International Mechanical Engineering Congress and 
Exposition under the auspicies of the American Society of Mechanical Engineers. The 
following sessions were organized and chaired by the P.I.: 

1. Heat Transfer in Microgravity Systems, 29th National Heat Transfer Conference, 
Atlanta, August 8-11, 1993. 

2. Heat Transfer in Microgravity Systems, International Mechanical Engineering Congress 
& Exposition, Chicago, November 6-11, 1994. 

3. Heat Transfer in Microgravity Systems, 30th National Heat Transfer Conference, 
Portland, Oregon, August 6-8, 1995. 

4. Heat Transfer in Microgravity Systems, International Mechanical Engineering Congress 
& Exposition, San Francisco, November 11-17, 1995. 
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5. Heat Transfer in Microgravity Systems, International Mechanical Engineering Congress 
& Exposition, Atlanta, November 17-22, 1996. 


4 Publications & Presentations 

4.1 Journal Papers 

1. Sadhal, S.S., Trinh, E.H. & Wagner, P., “Unsteady Spot Heating of a Drop in a 
Microgravity Environment,” Microgravity Science and Technology , in press (1997). 

2. Zhao, H., Sadhal, S.S., & Trinh. E.H., “Singular Perturbation Analysis of an Acous- 
tically Levitated Sphere: Flow About the Velocity Node,” submitted, Phys. Fluids 
(1997). 

4.2 Conference Proceedings 

1. Gopinath, A. & Sadhal, S.S., “Thermoacoustic Streaming Effects from a Sphere 
Subject to Time- Periodic Temperature Disturbances,” accepted. Tenth International 
Heat & Mass Transfer Conference, Brighton, U.K., August 1994. 

2. Trinh, E.H. & Sadhal, S.S., “Acoustic Streaming and Ultrasonic Processing of Low 
Melting Point Materials,” 1994 Int. Mech. Eng. Congress, Chicago, November 
6-11, 1994. In Heat Transfer in Microgravity Systems 1994 , ASME HTD-Vol. 
290, pp 43-52. 

3. E. H. Trinh & S. S. Sadhal, “Visualization of Internal Flows in Differentially Heated 
Drops,” Gordon Conference on Gravitational Effects in Physico-Chemical Systems, 
July 1995. 

4.3 Book 

Sadhal, S.S., Ayyaswamy, P.S., & Chung, J.N., Transport Phenomena with Drops and 

Bubbles , Springer- Verlag, 520 pages, 1997, ISBN 0-387-94678-0. This monograph in- 
cludes a chapter dealing with thermocapillary phenomena with drops and bubbles. 

4.4 Edited Works 

1. S.S. Sadhal, A. Gopinath, P.D. Jones, J. Seyed-Yagoobi, & K.A. Woodbury, Pro- 
ceedings of the the ASME Heat Transfer Division - Volume 1, American Society 
of Mechanical Engineers, No. HTD— 332, ISBN 0-7918-1519-6 (1996). 
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2. S.S. Sadhal, A. Gopinath, P.H. Oosthuisen & A. Hashemi, Proceedings of the 30th 
National Heat Transfer Conference - Volume 3, American Society of Mechanical 
Engineers, No. HTD-305, ISBN 0-7918-1704-0 (1995). 

3. S.S. Sadhal & A. Gopinath, Heat Transfer in Microgravity Systems 1994. American 
Society of Mechanical Engineers, No. HTD-290, ISBN 0-7918- 1408-4 (1994). 

4. S.S. Sadhal & A. Hashemi, Heat Transfer in Microgravity Systems 1993 , 

American Society of Mechanical Engineers, No. HTD— 235, ISBN 0-7918- 1148-4 
(1993). 

4.5 Invited Lecture 

"Fluid Flow and Thermal Transport with Acoustically Levitated Spheres,” presented at 
“Gotas, burbujas y pellculas (Drops Bubbles, & Films) ”, Escuela de Fisico-Quimica de 
Fluidos, Santander, Spain, September 9-13, 1996. 
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S. S. Sadhal, E. H. Trinh and P. Wagner 

Thermocapillary Flows in a Drop With 
Unsteady Spot Heating in a Microgravity 
Environment 


The unsteady localized spot heating of a liquid drop under 
zero-g conditions is examined theoretically. This pertains to 
space experiments to measure thermal properties of materials 
and the purpose here is to predict the thermal behavior of 
such systems. Spot heating can be achieved by a laser beam 
focused on a small region of the drop surface . The present 
theoretical model deals with situations of weak Marangoni 
flows , whereby the thermal transport is conduction domi- 
nated. The heat flow in the drop is treated as unsteady while 
the surrounding gaseous region is considered to be quasi- 
steady. The ensuing thermally driven flow is analyzed in the 
Stokes regime. 


1 Introduction 



Since the availability of experimental facilities in space, 
studies in low-gravity environments have gained consider- 
able attention^ Much of the interest in the last decade 
started with some pioneering papers of Subramanian and 
Weinberg [1, 2]. One of the principal advantages of a low^ 
gravity situation is that deep undercooling of liquids is 
possible because of fluid being undisturbed. In this context, 
it is of fundamental interest to acquire information about 
the thermophysical properties of the liquids in undercooled 
states. To this end, various low-g experiments have been 
proposed. Among them is the idea of spot heating a spher- 
ical sample of the liquid under test in zero, or low gravity 
(see fig. 1). This idea was proposed by Subramanian (see [3, 
4]) more than a decade ago but analytical studies have only 
been recently conducted. A laser beam can be used to heat 
a small fraction of the drop surface. The thermal field thus 
developed will lead to Marangoni convection and the obser- 
vation of the fluid motion can be used to infer the thermal 
properties of the sample. This possible use of method has 
also been investigated by Shen and Khodadadi [5, 6]. 
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Since the low gravity significantly weakens the natural 
convection phenomenon, Marangoni convection becomes 
isolated and it is the dominant mode of thermally driven 
motion. The study of Marangoni convection under micro- 
gravity condition has been the subject of considerable study 
in the last several years, particularly in relation to drops 
and bubbles. Among the major theoretical works are the 
investigations by Annamalai et al. [3], Barton and Subrama- 
nian [7-9], Kim and Subramanian [10, 11], Merritt and 
Subramanian [12, 13], Meyyappan et al. [14-17], and 
Shankar and Subramanian [18]. They have also addressed a 
number of problems involving compound drops. For spher- 
ically concentric systems, Bauer and Eidel [19] obtained 
three-dimensional solutions for a class of thermocapillary 
flows. They examined various configurations including the 
case of a drop within a drop with a rigid outer boundary 
and full continuity at the inner interface. In addition, they 
treated the situation of a rigid inner boundary with a free 
outer interface. The thermocapillary migration of a drop 
heated by radiation on one side has been analyzed in the 
steady state by Oliver and DeWitt [20]. They obtained an 
expression for the migration velocity as a function of a 
radiative intensity. Reviews of work on thermocapillary 
flows in microgravity are available in [21-23]. 


laser 

beam 



Fig. I. A schematic of a spherical drop being spot heated 
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Much of the earlier work has focused on drops and 
bubbles in a linear temperature profile and there has been 
little effort in the direction of spot heating. The present 
analvsis is directed towards the study of the thermal and the 
fluid dynamical effects of spot heating on a liquid drop. The 
Marangoni convection is considered to be weak enough so 
that the convective transport of heat may be neglected. 

Thus pure conduction dominates in the liquid phase, as 
well as in the surrounding gas phase. The drop is treated as 
fully transient, but the gas, owing to its low heat capacity, 
can be analyzed in the qua(sisteady state. This assumption 
limits the validity to cases for which the liquid thermal 
diffusivitv is much lower than that of the gas. For most 
glass-like materials in a gas, this is often the case. For 
example, the thermal diffusivity of glass is approximately 
6 10‘ 7 nvs" 1 , while that for air at 20 °C is about 
2 10’ 3 m 2 s~‘. The fluid flow is considered to be inertialess 
due to the high viscosity of the dispersed phase together 
with the weak thermal driving force. This corresponds to a 
low Marangoni number which is defined as 



where £) is the power of the laser; x, and fi correspond to 
the thermal conductivity, the thermal diffusivity, and the 
kinematic viscosity, respectively, of the liquid drop; and a is 
the interfacial tension. For glassmelts, we take the estimate 
given in table 1 of the thermophysical parameters near the 
solidification state [24]. 

Here, the upper limit of the viscosity is near the ‘transfor- 
mation temperature* in an undercooled state [24], Taking fi 
to be 10 3 N s m -2 , for Ma to be between the range 10~ 2 to 
10 -3 , £) should be about 0.6 to 6 mW. This power level is 
enough to set off an observable motion and yet keep 
Ma « 1. At higher viscosity values of the melt, the power 
level can be raised to obtain an observable motion. 

In the gaseous continuous phase, while the dynamic 
viscosity may be low, the kinematic viscosity (v = p/q) is 
large because of the low gas-phase density. The flow field is 
considered to be qua^s)eady under the low Reynolds num- 
ber approximation, i.e., the time derivatives in the Navier- 
Stokes equation are dropped. This is usually done for 
Stokes flow because the diffusion time scale, /? 2 /v, is much 
smaller than the dynamic time scale, R/U. This is consistent 
with the assumption. Re « 1. The drop deformation is also 
considered to be negligible because the normal stress forces 
arising from the surface tension variation and the fluid 
motion are similar in magnitude as the Marangoni forces. 

The flow observations coupled with some thermal mea- 
surements near the drop will provide sufficient data to 
acquire information about the thermal properties of the 
liquid under test. The transient evolution of the quasi- 
steady flow field prediction may be compared with the 
experimental measurements and thermal diffusivity infor- 


mation can be derived. The dominant transient is the heat 
diffusion in the liquid drop. This assumption realistically 
limits the validity to cases where i « x, v, v. The observa- 
tion of the evolution of the central ring of the toroidal 
vortex will provide information about the thermal diffusiv- 
ity. While the flow equations are tune-dependent, the tem- 
poral evolution appears in a quari-sjbady manner, as just 
mentioned, because the interface conditions are time-depen- 
dent. 


2 Statement of Problem 


The mathematical model consists of a liquid drop of radius 
R surrounded by an infinite gaseous medium. A laser beam 
is focused on a small area of the drop surface. The gas is 
considered to be transparent to the radiation which is 
absorbed at the surface with some reflected back and negli- 
gible penetration into the liquid. Thus the heated spot 
behaves as a surface source of heat. In the absence of any 
other driving force, the problem is symmetric about the axis 
of the beam. As stated earlier, we assume that the thermal 
diffusivity of the liquid is much smaller than that of the gas. 

We use the notation where the quantities corresponding 
to the drop are referred to with a hat ( A ). Thus, assuming 
i « x and neglecting convective transport, we have the 
temperature fields described by 

i ef d 2 f 2 ef i e / . _ a f\ ... 

a a/ dr 2 r dr r 2 sin (9) d8 \ ° ^ dd )' 

and 


. e 2 r 2 er 1 a / . an 

dr 2 * r dr * r 2 sin (6) dd V " ^ 26 / 
The continuity conditions at the interface are 


(3) 


ru = n.* 

and 

_ L too cos (0) 0 < 0 < 00. 

* Or erj,.* jo 0 O <0<K, 

where q 0 is the net heat flux deposited by the laser beam. 
This is distributed somewhat unevenly because of the sur- 
face curvature and this effect introduces the factor cos (#). 
However, it is^minor issue from a practical standpoint 
because the laser intensity is seldom uniform. 

In the far field, the temperature is uniform, i.e., 


(4) 

(5) 


r-r x as r x. 


( 6 ) 


The initial condition is 


f = 7 = 7 X at / = 0. (7) 

Tlfe differential equations for the fluid velocity can be 
reduced to a scalar form in the usual manner with the use 


Table /. Properties of liquid glass near the solidification state 



£ 

fi 

2 

*0 

da /d7 


l.OWm-' K- 1 

10 3 — 1 0' 23 N'sm" ! 

6^J0~ 7 m 2 s- 1 

300 N T m-’ 

-10- J Nm-' K- 1 


Microgravity sci. technol. X/l (1997) 


S. S. Sadhal et al.: Thermocapillary Flows in a Drop With Unsteady Spot Heating 


of the Stokes stream functions, i.e. f 
1 d\p -10^ 


U ' = r 1 
and 


r* sin (0) 00 * 


u r = • 


0l£ 


1 sin (0) 00 ’ 


U 9 r sin (0) dr * 

-1 0i£ 

= 


r sin (0) 0r ’ 
The stream functions, t p and \J/ obey 
Li, \p =0, and Li, I = 0, 
respectively. Here. 

5 2 sin (9) 5 ( 1 d 


L.,= 


dr 2 


59 \sin (9) 89 


( 8 ) 

(9) 

( 10 ) 

(H) 


is the Stokes operator. On the interface r = R, normal 
velocity is zero, and the tangential velocity is continuous. 
There is a tangential stress discontinuity that is equal to the 
surface tension gradient, ( \j R)(8al89). These interface con- 
ditions may be written as 


= fl-H = o. 


dip 

07 

and 


d\f 
: 07 


( 12 ) 

(13) 


Id 1 (dT\ , 

RdT\89j rmF ^ 

__ |~ r 0 / 1 d\p \ * r 

[ ^ sin (0) dr \r 2 dr J ^ sin (0) dr \r 2 dr 


a (\ 


d\f 


Jr-* 


(14) 

where p is the gas phase viscosity and p is that for the 
liquid. In addition, we must have finite velocity at the 
origin. This is satisfied by 


-“i ft < cc as r ■ 
r l 

3 Solution 


^0. 


(15) 


The solution consists of the determination of the tempera- 
ture field and the ensuing velocity field. 

J. 1 The Temperature Field 

The solution for the temperature distribution can be easily 
obtained by Laplace transform. Taking the Laplace trans- 
form of the set of differentia] equations (l)-(6) we obtain: 


&T 2 dT I_0 

0r 2 r dr + r 2 dr t 

d 2 T 2 dT 1 2 , . dT 

— 2 H i h ( 1 — *l 2 ) jr~ 

2 r dr r 2 dfj dr] 


\(sf- T^) = -rj 2 ) 

y. dr 2 r dr r i drj 


0 = 


dT 

dn 


dr 





L df 

r*e7 

an 


drj 



[ion 


n 0 <n < 1 
-i <n <no> 


(16) 

(17) 

(18) 

(19) 


where rj = cos (9), q 0 = cos(0) o and d>* = £//e. The solutions 
to these differential equations can be written as 


f — I a„i n (qr)PM + — , 

nm 0 S 


, £ (R Y +1 


i n (qR)PM+—> 

s 


( 20 ) 

( 21 ) 


where i n (qr) represents modified spherical Bessel functions 
and 


9 = 


1/2 


To satisfy the boundary condition (19) we expand the right 
hand side, so that 

(qR) + iSqR) " £ £ H * P M* 


I aA<t> k qi‘„ 

n - 0 l 


n -0 


where 


i 

= (n f i) J rjPM dn 

no 

_ (n -I- j )( 1 - rjo)\r± qPJjJq) - PM 0 )] 

(w*2X«-l) 

As a result we obtain 

n= 


, (n * 1), 


JssE\ 

V ks j 


J [q<P k Rr„(qR) + (n + l)i H (qR)] 
Therefore, we have 


. n) = ^ 


qpR \ 

ks ) m t o [q<j> k Ri‘ n (qR) + in * 1 )«„(**)] ’ 5 


HJJjr)PM 


fir 

and 

T, = (i2*\ v 1* 

(r ' n) \ ks ) .to [q<fi> Ri 'M R ) + (» + O'. (9 r )J s 

Inversion of the Laplace transform yields 

T - r --(ir) 


y [ y ^ k {j n {^r)l Jn (/., 

„.o U-I ( w 


,7?)]e 


-mL 


(rlRy 


n<j> k •*■(« + !) 




+ ])] - ;.l„R 2 4, 2 

PAcos (9)), (22) 


and 

9 


-r.-m 


y I y 2 4 > k e-*^-' 

1 ’ i (B - 1 )[nd>l + </>*-(« + 1)] - ;.L Rl <t>l 


n<P* 


i r ^}(f)"'«.Meos( 8 ),. (23) 
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where /. mw represents the charactenstic values given by the 
roots of 

- =("-*■ >)/,('■..»*). (24) 

and j n (x) represents spherical Bessel functions. 

3.2 The Velocity Field 

The stream functions represented by the solutions of equa- 
tion (10) are given as [25] 


[Mr. 6). fa, 9)] = X + 


find that 


PM dx. 


(25) 


co* (0) 


Upon satisfying the boundary and interface conditions 
(12), (13), and (15), the above general solution reduces to 


:or-(r j 


ftr. 6) = X CM 

n * l 

and 

' >m >- .?,«• {(*)““ -(s 


/>„(.x) dx (26) 


CO* (ff) 


PM) dx. (27) 


CO* (0) 


Applying the tangential stress jump condition, 

1 S<7 

Rd8 = ^ r< ~ 

_ _r r a / 1 df\ . r 5 / 1 Stfr 

L" sin (0) dr \r 2 dr) M sin (6) Sr \r 2 Sr 


- * 

(28) 


we obtain 
da 2(n (3) 


S6 /? 

Now, using 

da _ dff /an 

66 _ drUeA.* 

and 


, f C„(0(2n h- 1) - ;‘ n p;(cos (fl)). (29) 


■+■ 1) 


a., 


(30) 


3 T\ (2 qoR\ 

ce) r . K l k J 

« X 

• z z 


/y.^(cos(fl))(-sin(fl))^e- M i- t 
I „tr 0 (« -*-4>k -(« + 1)] 

(31) 


i /y,f:(cos(fl)X-Hn(g )) 

£ « - o n<t> k ■+■ (n 1) 


r i 


|_ 2 0 < - 

[dr) 

f * 

J V 

2e-“~ 


U-i (« * 1)[« +</>*- <££(« - 1 - 1)] - 

1 1„ [ »(«•*• 1) 1 
rt *+■ <£*(* l)j "L 2* h- 1 J 

J.J Thermocapillary Force and Drop Migration 


( 32 ) 


The force experienced by the drop is proportional to the 
strength of the stokeslet, and is given by 


F = 


4rt^C, 


which has the steady state value 




2nq 0 R 2 ( 1 — ^p) 

3(1^)(2 + <M*J 


(-£> 


(33) 


(34) 


Therefore, the thermocapillary force would cause the drop 
to move, unless an additional neutralizing force is present. 
Such a force can come about by means of a weak acoustic 
field, or an electrostatic field. For such a situation it is, of 
course, understood that the secondary flows arising from 
the neutralizing force are negligible, and the consequent 
limits of validity apply. In the absence of such a force, the 
migration velocity can be obtained by simply adding the 
Rybczynski-Hadamard velocity field to the solution given 
by eqs. (26) and (27). That is, for a migration velocity, L\ 
we add on 


*o(r. 0) = \ UP 2 


2(1 


((;)'-© 

«[(*)“( ')]} !in,(< ’ ) - 


2 — 30 


(35) 


and 


8 > - J - (s)'] 5 ™’ <*>• 1361 

to i£(r, 6) and \J/(r , 0). respectively, and require that the net 
force on the system be equal to zero. This is achieved by 
letting the term proportional to (r sin 2 (0)) (planish. As a 
result, the following expression for the migration velocity is 
obtained: 


, C x [ 2(1 
' * R 2 l 2 + 3<j> M _ ' 


2^3 

It has the steady state value, 

U„ = 


q a R(\-r]l) 


_lkn{2 3<£„)(2 

{ dr) 


(37) 


(38) 


which is consistent with the result of Oliver and DeWitt [20] 
for the special case of 0 O = 


4 
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3.4 Drop Deformation 

Upon integrating eq. (29), the surface tension distribution 
is found to be 


/m 20 1 + P) £ ^ 

ff(9) =<J 0 — 2 — 2- c *( l > 


' 2n + I 
n(n + 1) 


P n (cos(6)y (39) 


Following Sadhal , Ayyaswamy and Chung [21], the defor- 
mation may be expressed as 

R(0) = Ro[ l+CiO)]. (40) 

where 

Am L 


r(Q\ = 2(/i + £) “ CO) 

„. 2 n(n 1) 


.(«+2) («~2). 


/»„(cos(0)). 


(41) 

and /? 0 is the radius of the isothermal drop. The surface 
tension of the liquid at the reference isothermal state is 
denoted by < 7 0 . By using the expression for C„(/), it is not 
difficult to see that 


rta\ — 9 ° 

a6) -^j[df) 

oo C ac 

■ I I 7 

<i ■ 2 Ow - 1 v 


2c 




1 


1)[« + 1)1-2 
1 




>K 


n + <f> k (n + \)) * L(n + 2) (/! - 2)J 


^(cos(fl)). 
(42) 

Here, considering £=0.1 to lOmW, a beam radius of 
1 mm, a drop radius R 0 = 5 mm, and other properties as 
given in table 1, the dimensionless parameter in front of the 
above expression for ~{9) is approximately given by 

q 0 R 0 /dff\ 


<x 0 £ \dTj 


U10-’ to IOC 


(43) 


It is therefore easy to see that for the estimated set of glass 
properties, the drop deformation is very small. 


4 Results and Discussion 


The important results from this analysis consist of the 
temperature field and the flow streamlines. A computer 
program was written in FORTRAN to run on the Cray 
Y-MP at JPL. The program carries out the formulation in 
stages. At the first stage the roots, A mn , of the characteristic 
eq. (24) are obtained by the regula falsi method. For 
accuracy, the Legendre polynomials and their derivatives 
are evaluated using recurrence relations at higher orders 
rather than using polynomial expressions. The spherical 
Bessel functions are treated differently according to 
whether 

( 1) the square of the argument is very small compared with 
the order, in which case an expansion in this ratio is 
used; 

(2) the order is greater than 3, in which case backward 
recurrence from an order based on the desired accuracy 
is used; 

(3) otherwise, simple expressions with sines and cosines are 
used. 


Because more terms are required for convergence of the 
steadi-state solution than for convergence of the transient 
response, each is tested separately. The results are displayed 
as isotherms and streamlines using standard NCAR graph- 
ics library routines. 

The temperature is scaled with the parameter, q^Rjlc, and 
the calculations have been carried out at Fourier number 
values, Fo = itjR 2 = 0.1, 0.2, 0.5, 1, 2, 10. The heated spot 
is defined by the angle, 0 O = 15°, as show-n in fig. 1. For 
each case, the isotherms and the flow streamlines are given 
in figs. 2a-f. The diffusion of the heat can be seen progress- 
ing in time, as one set of isotherms is compared with the 
next. 

The flow streamlines exhibit a single vortex arising purely 
from the thermally driven motion. The center of the vortex 
is away from the equatorial plane and finds a position 
toward the heated spot. A similar shift was found by Oliver 



Af g. 2. Flow streamlines for a spot-heated drop, <p k = 100, 6 0 - IS 1 , 
(a) Fo = 0. 1, (b) Fo = 0.2, (c) Fo = 0.50 , (d) Fo = 1.0, {e) Fo = 2.0, 
c/5 Fo = 1 0.0 
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and DeWiu [20] for the case of a drop heated by radiation 
over half of the free surface. In the present analysis we find 
that as the thermal diffusion progresses, the vortex center 
shifts towards the equator. The position of the vortex 
center can be used to infer the thermal diffusivity of the test 
liquid. The usefulness of this approach is particularly evi- 
dent because the streamlines are independent of the liquid 
and the gas viscosities at the present level of approximation. 



List of Symbols 

(„(*) modified spherical Bessel function 

_/„(*) spherical Bessel function 

L_, Stokes operator 

Ma Marangoni number, Q( -da IdT)l(kpi) 

P„(x) Legendre polynomials 

Q laser beam power [W] 

q (sji)" 2 

q Q heat flux 

r radial coordinate 

R radius of the drop 

s Laplace transform parameter 

t time 

T temperature ; 

u r , u g velocity components 

2 thermal diffusivity 

6 polar angle 

0 o angle subtended by heated area 

\x dynamic viscosity 

*7 cos (0) 

o interfacial tension 

r* shear stress I 

<t> k conductivity ratio, (fcfk) 

\j/ stream function 

x far field (subscript) 

T Laplace transform of T 

( A ) dipersed phase quantities ^ 

1 


Acknowledgements 

Support of this work by the U.S. National Aeronautics and Space 
Administration (Grant Nos. NAGW-3378 and NAG3-1842) is 
gratefully acknowledged. 

References 

1 Subramanian , R S , Weinberg, M C. : The Role of Convective 
Transport in Dissolution or Growth of a Gas Bubble. J. Chem. 
Phys., vol. 72, p. 6811 (1980) 

2 Subramanian , R S., Weinberg, M C. Asymptotic Expansions 
for the Description of Gas Bubble Dissolution and Growth. 
AIChE Journal, vol. 27, p. 739 (1981) 

3 Annamalai, P ., Shankar , N., Cole , R., Subramanian, R. S.: 
Bubble Migration Inside a Liquid Drop in a Space Laboratory. 
Appl. Sci. Res., vol. 38, p. 179 

4 Subramanian, R. S.: Applications in the Academic and Scien- 

tific Community. IBM: Perspectives in Computing, Vol. 4, No. 
2/3 (1984) w - 


5 Shen, F., Khodadadi , J. M.: Transient Marangoni Convection 
Inside a Liquid Droplet. APS Bulletin, vol. 39, p. 1841 (1994) 

6 Shen, F., Khodadadi, J. M.: An Extended Laser Flash Tech- 
nique for Thermal Diffusivity Measurement of High Tempera- 
ture Materials; in: NASA Conference Publication 10121, 1992. 
Workshop on Thermophysical Properties of Molten Materials, 
Cleveland, OH, p. 203 (1993) 

7 Barton, K. D., Subramanian, R. S.: Thermocapillary Migration 
of a Liquid Drop Normal to a Plane Surface. J. Colloid 
Interface Sci., vol. 137, p. 170 (1990) 

8 Barton, K. D., Subramanian, R. S.: The Migration of Liquid 
Drops in a Vertical Temperature Gradient. J. Colloid Interface 
Sci., vol. 133, p. 211 (1989) 

9 Barton, K. D., Subramanian, R. S.: Migration of Liquid Drops 
in a Vertical Temperature Gradient - Intersection Effects Near 
a Horizontal Surface. J. Colloid Interface Sci., vol. 141, p, 146 
(1991) 

10 Kim , H. S., Subramanian , R. S.: Thermocapillary Migration of 
a Droplet With Insoluble Surfactant. I: Surfactant Cap. J. 
Colloid Interface Sci,, vol. 127, p. 417 (1989) 

1 1 Kim, H. S., Subramanian, R. S.: Thermocapillary Migration of 
a Droplet With Insoluble Surfactant. II: General Case. J. 
Colloid Interface Sci., vol. 130, p. 112 (1989) 

12 Merritt, R. M., Subramanian, R . S’.: The Migration of Isolated 
Gas Bubbles in a Vertical Temperature Gradient. J. Colloid 
Interface Sci., vol. 125, p. 333 (1988) 

13 Merritt, R. M., Subramanian, R. S.: Migration of a Gas Bubble 
Normal to a Plane Horizontal Surface In a Vertical Tempera- 
ture Gradient. J. Colloid Interface Sci., vol. 131, p. 514 (1989) 

14 Meyyappan, M , Subramanian , R. S Thermocapillary Migra- 
tion of a Gas Bubble In an Arbitrary Direction With Respect 
to a Plane Surface. J. Colloid Interface Sci., vol. 1 15(1), p. 206 
(1987) 

15 Meyyappan, M., Subramanian , R . 5.; The Thermocapillary 
Motion of Two Bubbles Oriented Arbitrarily Relative to a 
Thermal Gradient, J. Colloid Interface Sci., vol. 97(1), p. 291 
(1984) 

16 Meyyappan , M., Wilcox, W, R., Subramanian, R. S’.. The Slow 
Axisymmetric Motion of Two Bubbles in a Thermal Gradient. 
J. Colloid Interface Sci., vol. 94(1), p. 243 (1983) 

17 Meyyappan , M., Wilcox, W. R., Subramanian , R. S.: Thermo- 
capillary Migration of a Bubble Normal to a Plane Surface J, 
Colloid Interface Sci., vol. 83, p. 199 (1981) 

18 Shankar, M., Subramanian , R. S.: The Slow Axisymmetric 
Thermocapillary Migration of an Eccentrically Placed Bubble 
Tnside a Drop in Zero Gravity. J. Colloid Interface Sci., vol. 
94(1), p. 258 (1983) 

19 Bauer, H. F., Eidel , W,; Marangoni-convection in a Spherical 
Liquid System. Acta Astronautica, vol. 15, p. 275 (1987) 

20 Oliver, D. L. R., DeWitt, K. J Surface Tension Driven Flows 
in a Micro-gravity Environment. Int. J. Heat Transfer, vol 31. 
p. 1534 (1988) 

21 Sadhal, S. S., Avyaswamy, P. S., Chung, J. N.: Transport 
Phenomena with Drops and Bubbles. Spnnger-Verlag. New' 
York (1997) 

22 Subramanian, R. S.: The Motion of Bubbles and Drops in 
Reduced Gravity; in: Transport Process with Drops and Bub- 
bles, R . P. Chhabra , D. DeKee (eds.) p 1-32, Hemisphere 
(1982) 

23 Wozniak, G., Siekmann, J., Srulijes, J.: Thermocapillary Bub- 
ble and Drop Dynamics Under Reduced Gravity - Survey and 
Prospects. Z. Flugwiss. Weltraumforsch, vol. 12, p. 137 (1988) 

24 Schohe, H.: Glass. Nature, Structure and Properties (trans- 
lated by M. J. Lakin). Spnnger-Verlag, New York (1990) 

25 Happel , J ,, Brenner, H.: Low Reynolds Number Hydrodynam- 
ics. Martinus Nijhoff, The Hague (1983) 



6 


Microgravity sci. technol. X/l (1997) 



Singular Perturbation Analysis of an Acoustically 
Levitated Sphere: Flow About the Velocity Node 
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Abstract: This analysis consists of the development of the fluid flow about a spherical 
particle placed at the velocity node of a standing wave. High-frequency acoustic fields are 
being used to levitate particles in Earth gravity, and to stabilize particles in low-gravity 
situations. While a standing wave in an infinite medium may be purely oscillatory with 
no net flow components, the interaction with particles or solid walls leads to nonlinear 
efffects that create a net steady component of the flow. In the present development, the 
perturbation procedure of Riley [9] is employed to derive the flow field for the situation 
when a spherical particle is positioned at the velocity node. As in Riley s [9] solution 
which applies to a sphere at the velocity antinode, it is found that there is a thin shear- 
wave region adjacent to the spherical boundary. However, this thin Stokes layer does 
not cover the entire sphere as in the previous case, but lies mostly around the equatorial 
region of the sphere. On the remaining surface of the sphere the Stokes layer opens to 
the surrounding flow field. 
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Nomenclature 


a radius of the sphere 
A/ 2 frequency parameter, n 2 jj/v 
r radial coordinate, r* / a 
R Reynolds number, U 00 n/i' 

R s streaming Reynolds number, U 2 JuJV = £ 2 M 2 
t time 

Uoo velocity amplitude in the sound wave 
2 coordinate along the polar axis 


Greek Letters 

6 Stokes layer thickness, \jv/u) 
e perturbation parameter, £ = U^jua 
t) normal boundary layer coordinate [equation (18)] 
9 angular coordinate 

A wavelength of the sound field, viscosity ratio 

[l cos 9 

v kinematic viscosity of the fluid 
q radial cylindrical coordinate, q = rsinO 
r dimensionless time, u it 
i>, 'I' dimensionless stream function 
uj angular frequency 


Superscripts 

* represents nondimensional quantities 
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1 Introduction 


The potential for the development of both fundamental and applied sciences in the area 
of liquid undercooling with containerless experiments is now well established. With 
the achievement of deep undercooling in zero gravity, the possibility of new materials 
is becoming evident. It is now well known that this high degree of undercooling can 
promote certain types of crystal growth and at the same type provide homogeneity of the 
product. 

In an Earth-gravity field, some degree of undercooling can be achieved by acoustic 
and/or electrostatic levitation. With the application of the principle of radiation pressure, 
ultrasound levitators have been in use for many years in ground based experiments [1, 
2, 3, 7], The recent advances in single particle levitation technology have rekindled 
scientific interest in the determination of the bulk and surface physical properties of 
liquids in the metastable supercooled, supersaturated, or superheated states. The capability 
of indefinitely suspending in a host fluid or vacuum, a high-purity liquid in the form of 
a free drop, without inducing large scale translational or oscillatory motion has created 
the opportunity to accurately determine the temperature dependence of a number of its 
thermodynamic properties. It has also become possible to observe and to quantify the rates 
of transport processes involving its free surface and the surrounding fluids, and to examine 
the details of phase transformations such as melting and solidification, evaporation and 
condensation. 

The acoustic field provides the radiation pressure necessary to levitate a liquid drop in 
a gravitational field. The studies on the effects of radiation pressure on spheres and disks 
goes as far back as the 1 930 when some of the earliest theoretical studies were conducted 
by King [4, 5]. With the application of this principle, ultrasound levitators have been in 
use for many years in ground based experiments. Electrostatic levitation is another way 
for containerless processing. In zero gravity, the acoustic field can be used to stabilize a 
fluid particle. 

Acoustic streaming can be classified as two common types. One happens because of 
the spatial attenuation of a wave in free space, e.g., an attenuating beam of plane traveling 
wave. This type of streaming is usually associated with a high Reynolds number flow. 
The second mechanism arises from the friction between the medium and a solid wall 
when the former is vibrating in contact with the latter, e.g., a wave traveling down a 
wave-guide, a standing wave in a resonant chamber, or a wave scattering off a solid 
object. Unlike the spatial attenuation mentioned earlier, this effect is largely confined 
to a thin viscous boundary layer of thickness 6 = (2u/u>) 1 ^ 2 on the surface, where u is 
the kinematic viscosity of the medium and u is the angular frequency of the wave. It is 
also a significant dissipation mechanism, and provides a strong force in driving acoustic 
streaming. While the medium outside the layer vibrates irrotationally as in a sound field, 
the one within the layer has vorticity because its motion has to conform to the no slip 
condition on wall. We are now interested in this second mechanism of streaming. 

If a body of typical dimension a oscillates with velocity, say, Uoo cos (u>t) in a viscous 
fluid and e — U^juja -C 1, then, although the leading order solution is oscillatory, higher 
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node 



antinode 


Figure i: Velocity node and antinode identification 


order terms include not only higher harmonics but steady contributions to the velocity. 
Mathematically, this can be explained by existence of the nonlinear terms which may have 
steady nonzero component. Physically, the condition e<l implies that the amplitude of 
the oscillation is small compared with the particle radius a. The existence of this steady 
streaming was first pointed out by Rayleigh [8] in his work on Kundt’s dust tube and was 
later studied in a boundary layer context by Schlichting [10] who considered flows with 
the additional constraint |M| 2 = wa 2 / v 1, where u denotes the kinematic viscosity 
of the fluid. For such a flow it is now well established that the first order fluctuation 
vorticity is confined to a shear- wave region of thickness 0(v/lj) 1/ 2 beyond which steady 
velocities 0{eU 00 ) persist. Riley [9] has considered the case of an oscillating sphere for 
both \M\ » 1 and 0 < \M\ < 1. He calculated the streaming around the sphere which is 
at the velocity antinode of the wave that vibrates vertically. 

In the analytical investigations to date, little consideration has been given to a particle 
placed at the velocity node of the wave. The velocity node and antinode are identified 
in Figure 1 . The results from the investigation of the flow field will provide informa- 
tion about the characteristics of the levitation process. The theory will be very useful in 
overcoming some of experimental problems by providing suitable direction. For exam- 
ple, presently with acoustic levitation there is a residual flow field, including solid-body 
rotation for drops. Such instabilities interfere with, for example, the steady application 
of thermal stimuli which would serve as a means for the measurement of thermophysi- 
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cal properties. This problem needs to be solved and detailed understanding of the flow 
studies would be beneficial. Furthermore, very little is known about the pressure field. 
The flow calculations together with experimental measurements of the drop position can 
be used to predict the pressure field which is difficult to measure. For levitation under 
zero-gravity conditions, the drop assumes an equilibrium position at the velocity antinode 
when the external medium is a gas. When the particle phase has higher compressibility 
than the external phase (e.g., a gas bubble in a liquid), of if the particle is very small, 
the equilibrium position is at the velocity node. While the antinode solution has been 
available from Riley’s [9] work, the node solution is a new development. 

In the discussion here, the focus is on the analysis of solid sphere being placed at 
the velocity node of the wave. Besides the direct application to such physical cases, the 
analysis about the velocity node leads to an important result for calculating the streaming 
when the sphere is placed between the velocity node and antinode of the wave. Eventually, 
the analysis will be used to obtain the flow field for a sphere levitated at a position 
between the node and the antinode. In this proposed development, use will be made 
of the existing antinode solution of Riley [9] and the present node solution through a 
nonlinear combination. Although Lee & Wang [6] have considered this kind of problem, 
their result depends on an algorithm for calculating the tangential velocity on the edge 
of the recirculating shear layer. Detailed flow field in the shear layer for such situations 
are yet unavailable. A more rigorous development of the flow field is necessary to 
fully understand the fluid mechanics, and the present calculations provide an additional 
component for the complete development. 

Riley [9] gave the solution for an oscillating sphere in an otherwise quiescent infinite 
fluid medium. This solution is applicable to to a small sphere positioned at the velocity 
antinode of a standing wave provided a«A, i.e., the size of the sphere is small compared 
to the wavelength. We briefly discuss Riley’s [9] solution since it forms a basis for further 
development. The following dimensionless parameters are relevant: 



v 


M 2 


iua 2 , 

and 

v 


R 

\ Mp 


Uoo 

= £ = <sC 1, 

ua 


( 1 ) 


where is a characteristic velocity and R is the Reynolds number. While Riley [9] 
considered both \M\ 1 and \M\ » 1, the latter case (high frequency) is the one relevant 

to ultrasound levitation. 

For a standing wave with velocity 


u z = U'^ sin kz e lwt , 


( 2 ) 


the local velocity in the neighborhood of the node (z = 0) is 

Uz = U'oo (^ z ~ ^ el0Jt • ( 3 > 

With a small particle at the antinode, the surrounding field may just be taken as the first 
term 

u 2 =0^ ,Wt = ^Wa)e w ‘ (4) 
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where — U'^ka. 


2 Equations of Motion 


By scaling the flow parameters as follows: 

u = ip* = jr— x* = -, and r = ut, (5) 

b OC V ryo( l (l 

and dropping the asterisks, the Navier-Stokes equation of motion, in the stream function 
formulation, may be written as 


d(D^P) £ 

dr r 2 


d(ip, Dhp) 
d(rji) 


+ 2 D 2 ip Lip 


D A ip 

]mF’ 


( 6 ) 


where 

»■ - 

and fi - cos 0. The Reynolds number R, the frequency parameter M, and the perturbation 
parameter e are defined in equation (1). 

The stream function ip is related to the velocity components as follows: 


u r 


1 dip 
r 2 dfi 


and 


U g = — 


(1 -fi 2 ) 5 dip 
r dr 


(7) 


With this formulation the continuity condition for an incompressible fluid is satisfied. The 
flow may be regarded as incompressible for low Mach number. 

The boundary conditions are 


iP = 


dip 

Ik 


= 0 on r = 1 


(no slip), 


( 8 ) 


and 

ip — ► kg 2 ze iut = I r 3 /z(l - p?)e rT as r -* oo, (9) 

where g = r sin 9. This far-field condition is derived from the velocity field in the 
neighborhood of the node, given by equation (4). 
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3 Solution 


We seek a perturbation solution in the in the outer flow region by writing 

w = 0o + £0i + 0(s 2 ) (10) 

Similarly, in the shear-wave layer the following perturbation expansion is used: 

^ = vkg T £ 4* i -1 ~ O (s 2 ) , (11) 


where 'I' is related to 0 through 



( 12 ) 


Here R s = sR is the streaming Reynolds number which is also considered to be small. 
This relationship (12) indicates that the appropriate scale for the stream function in the 
shear- wave layer is U 00 o.(u/u>)^ rather that U^a 2 . 


3.1 The Leading-Order Solution 


Upon inserting equation (10) into equation (6), we find that the leading order stream 
function for the outer flow satisfies 


I M | 2 


d(D' l ip q) 
Or 


£> 4 0 0- 


( 13 ) 


For | M | 1, the vorticity, generated at the sphere, is confined to a thin ‘shear- wave’ 

layer of thickness 0 ( | M | ~ 1 ) . Outside this thin layer the flow field is irrotational and 
described by 


djDhM 

dr 


( 14 ) 


The boundary conditions are 


0o = 0 on r = 1 


and 0o ~ |A:ar 3 /i(l — fi 2 )e tT . 


( 15 ) 


The solution to equation (14) with the same behavior in /x as the far field is found to be 


0o = \ka(r 3 - r 2 )/i( 1 - /x 2 )e ir . 


( 16 ) 


For the inner region, we substitute the expansion (11) into the Navier-Stokes equation 
(6) and obtain for the leading order. 


where 


d 3 *o = 1 /d 4 fr 0 \ 

drdrj 2 2 \ dij 4 ) ’ 


n = (r - 1) 


|M| 

V2' 


( 17 ) 

( 18 ) 
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The boundary conditions for * f) are 


dr 7 


= *o = 0 


at 


V = o. 


The solution is found to be 


(19) 


$0 = c {(1 + i)ri - 1 + e- (1+,:) ’'} fi( 1 - fi 2 )e' T , (20) 

where c is an undetermined constant. Upon the expansion of do about r = 1, and matching 

#o with '0o when rj — > oo, we obtain c = |A;a(l — i). 

3.2 The First-Order Solutions [0(f)] 

The first order solutions are much more complex than the leading order solutions. With 
the substitution of the expansion given by equation (11) into the momentum equation (6), 
we obtain the following differential equation for * ( : 

<9* 0 / \ <9*0 /<9 3 * 0 \ / 2/“i \ d 2 * 0 /<9 *o\ _ 1 /<9 4 *i\ 

drdr] 2 drj \dpdr] 2 J dp \ dr / 3 j yl — ft 2 ) dr] 2 \ dr] ) 2 y dr / 4 j 

( 21 ) 

The solution to this equation may be written as a two-term eigenfunction expansion, 

*! = (ka) 2 [/i(?7, r)/2(l - /l 2 ) + / 2 (t?, r)/l 3 (l - /i 2 )] . (22) 

By the substitution of (22) into (21), we find 


1 

2 


'd*f i' 
dr? 4 


]i{\ -n 2 ) 


and 

1 (&h\ 

2 \ dr ] 4 J 


/i 3 (l-/i 2 ) 


^2/1(1 - p 2 ) + f| [2ir]e-^e 2iT - 2rr]e~^e- 2 ' T 

+ 2iT]e- {1+i)r > + 2e- (1+i)r, (l - e - (1 ~ i),? ) - 2 
+2e- (1 - t)rj (l - e- (1+l)TJ )] p( 1 - p 2 ) (23) 

^ 3 (1 - p 2 ) - || [2ir]e~^ i)ri e 2lT - 2ir]e~ {1 - l)ri e- 2iT 

+ 2ir } e~ (l+i)r ' + 2e~ {1 +^(l - _ 2 

+ 2e _ ^ _ * , ’ ? (l - cHi-wiu)} ^3(i _ p2^ 

+ f {(1 + i)e~ (1+i),? (l - e~( l+i)TI )e 2iT 
+ (1 - i)e-( 1 " < 0(l - e-V- i)T, )e~ 2 ' T 
+ (1 — i)e _ 0 _i),, (l — e ~( 1+i 0) 

+ (1 + i)e- (1+l) "(l - e-0- ),? )} p 3 ( 1 - p 2 ). (24) 
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By separating the steady streaming from the unsteady part, we can write fi as 

fi{V,T) = Qoin) + Qi 2(v)e 2lT - 

With the insertion of (25) into (23), it is found that 


£ / 

2 \ dif 


f + e" (1+ ' )?? (l - e- (1 “ i)7 ') - ir}e~ {x - i)ri 


+e 


-(1 ~i)Vi 


1 — e 




,. )V) ? 


and 


j = i~f + fine~" +i) \ 

4 \ dr) 4 j dr} 2 8 

with boundary conditions 

C20 = ^ = 0 and C22 = = 0 at r, = 0. 

dr} dr] 

The solution for equations (26) and (27) are found to be 

C20 = - |f e _2T/ - ~ r}e~ v sin r} - Y e ~ V cos *7 _ si» ^ + If ~ it 7 ? 

C22 = §§[\/2(l + i)je-^ 1+,) ’- |f[v/2(l+i)] + five- {1+l)v . 

As in the case of fi(rj,T), the function / 2 (r?, r) may be decomposed in the form 

h{r},r) = C 30 (7) + (32(v)e 2tT - 

By the substitution of (30) into (24), we obtain 


and 


d 4 Qo 

dr } 4 

1 ( ^ 4 C32 X 


Y —3 rje v sinr7 + 3e 2r> — 2e n cos 7} + e v sin rj 


= 1- 


. d 2 C,3i 


4 \ dr f ) dr f 

with boundary conditions 

and 


- + |(i + i) e -< 1+ *>^(l - c-t 1 ^), 


C 30 = = 0 at rj = 0 


C32 = 


dr} 

( K:u 

dr} 


0 at 77 = 0. 


( 25 ) 


(26) 


(27) 


(28) 

(29) 

(30) 

(31) 


(32) 
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The solutions to equations (31) and (32) are found to be 

Cio = w e ~ 2V + ? r ? e_ '' V + 25e - ' 7 cos q + sin q - ^ + j§q (33) 

<32 = f(l + i)e- (I+l) ' , -fi7?e- (1+7) ' 7 + ^(l + /)c-' 2 ( 1+ ' : >'' 

0 0 64 

+ ||(1 + i) [(llv^ - 9) - ll\/2e -v/2(1+ 9 f) j . (34) 

As a result, 

'I'j = (kn) 2 {(— §§e 2rJ — ^-qe~ v s'mq — ^e _,? cos77 — ^e - ' 7 sin q + ^ — y|j^ 

+ [|(1 + i)>/2 - l) + f iqe~^} e 2 "} ,1(1 - /l 2 ) 

+ { ( 5§ e_2T? + ^-qe~ n sinq + 25e~ n cos q + ^pe -77 sin q — ^ T 

+ [f (1 + - 2 iqe-^ + g(l + i)e ~ 2 < 1+ ^ 

+ §f (1 + i) (11^2 - 9e-^ (1+7),J )] e 2ir },l 3 ( 1 - ft 2 ). (35) 

Thus, the first order velocity has a steady part that exists outside the shear-wave layer 
and induces a steady streaming velocity O(s) outside that region. Upon letting q — ► oo, 
we see that at the edge of the shear-wave region, 

'P ~ |(fca) {t?cost — (|)\/2cos(t — ^7r)j ,1(1 — ,1 2 ) 

+ e( ka ) 2 {(^ - ft 7 ?) “ it cos ( 2r + H} £(1 - V 2 ) 

+ e i ka ? {(-§ + 'j§v) + |(ll\/2-9)cos(2r-|- f?r)}/i 3 (l -/l 2 ). (36) 
Here it is clear that 'P has the steady component 

= £(|)[(l3 - 35/7 2 )-6(l-3,- l 2 )^]A(l-/'< 2 ) + 0(£ 2 ). (37) 

In the outer region where r — 1 = 0(1), we need an asymptotic solution that should 
satisfy the far-field condition in equation (15), and as r — > 1, it should match with the 
inner solution as q — + oo. Therefore, in the outer region we have equation (10) for c 
where ipi satisfies the condition that as r — ► oc, 01 = o(r 3 ). The boundary condition 
in the inner region can be obtained from the matching condition and we may infer the 
behavior of in the inner region in terms of the inner variables [equations (18) and (12)] 
by letting r — * 1 in equation (36) for 'P. The equation for the function ip\ is obtained by 
equating coefficients of powers of e in (6). The terms of O(e) give 



(38) 

We can further write tp x as 


Vh = F } (r, ,1) + G, (r, /l)0(r), 

(39) 
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thus separating the steady streaming from the unsteady part. Upon inserting (39) into 
(38), we obtain the following equation for G x : 

D 2 G X = 0, (40) 


with boundary conditions 

G\ = o(r 3 ) as r — > oo 
Gi ~ p{\ — p 2 ) as r — + 1. 
The solution of (40) for Gi(r, p) is 


- 1 1 ' 2 ) 
r 2 


(41) 


As discussed above, the matching conditions for (36) and (39) determine <p(r). Therefore, 
letting T) — ► 0 in equation (36) and matching O(e) terms with it is not difficult to see 
that l 

4 >{ t ) = — 1( ka)Rs 2 cos(r — ^7r). 

Following Riley’s [9] perturbation method and using some scaling arguments, F x (r, /l) 
should satisfy 


D 4 F x (rJi) = 0. 


(42) 


To solve (42), an intermediate function H(r,n) is defined so that 

D 2 F x = H 
and 

D 2 H = 0. (43) 


The solution for H(r, /i) is 

00 _i 

H(rji) = Z(A n r” +1 + Bnr-nCj^) 

n = 0 


where 



(1 - fi 2 ) dPn 
n(l + n) dp, 


for n > 1 


(44) 


and 

c;Hp) = -p 

in which P n [p) is the Legendre polynomial. The solution for F x (r,p) takes the form 

00 _i 

F i( r -^) = 'j£ l fn(r)C n j i(/i) (45) 

n=0 
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where f n ( r ) satisfies 


fn(r) + \(-n(n+l))f n (r) = A n r n+l + B n r~ n . (46) 

r z 

Thus the solution to equation (46) is 

fn(r) = c n r n+l + D n r~ n + E n r n+ 3 + F n r~ n+2 (47) 

where A n , B n , C n , D n , E n and F n are constants to be determined. The function /„(r) 
must satisfy the condition that f n (r) = o{r 2 ) asr->oc and the inner boundary matching 
condition. So, F, (r, /l) can be written as 

F, (r, /l) = (D>r ~ 2 + F 2 )fi( 1 - /l 2 ) + (ZV ~ 4 + F 4 r~ 2 ) ( 7/1 3 - 3/7) ( 1 - /l 2 ) 

= ( D 2 r~ 2 4- F 2 — 3 F> 4 r -4 - 3F 4 r~ 2 )/7(l - /l. 2 ) 

+ 7(F> 4 r~ 4 + F 4 r- 2 )/l 3 ( 1 - £ 2 ). (48) 


Now by letting r — ♦ 1 and introducing the inner variable r/ defined in (18), we have 


F,(r, fl) ~ | (F> 2 + F 2 - 3F» 4 - 3F 4 ) + (-2 F> 2 T 12F> 4 + 6 F 4 )^r/ 1 /l(l - /l 2 ) 

l Fl J 


+ I (7£> 4 + 7F 4 ) + (—28 F> 4 - 14F 4 )^?r7 1 /1 3 (1 - ,7 2 ) (49) 

l FI J 

and matching to this order with (36), 

F>2 + F 2 - 3F> 4 - 3F 4 = 0 (50) 

-2D 2 + 12F» 4 + 6F 4 = -{§ (51) 

7£> 4 T 7F 4 = 0 (52) 

— 28£>4 — 14F 4 = ^>. (53) 


From equations (50) — (53), we deduce that 


D 2 


75 
1 12 1 


75 

112 ’ 


Da 


225 

224 ’ 


and F 4 = 


225 
224 * 


With these constants determined, the expression for i/’i is found to be 

•01 = -|(A:a)Fr 5 r _2 cos(r - - /i 2 ) 

+ ^(fca) 2 (l -M 2 ) 

+ §§ (fca) 2 (r -2 - r~ 4 )(7/i 3 - 3/i)(l - fi 2 ). (54) 


12 







4 Discussion 


The inner and outer solutions both have time-independent components which exhibit the 
steady streaming motion. The flow field exhibits symmetry about the polar axis and 
the equatorial plane. This symmetry is to be expected since the problem, as defined, is 
symmetric. The streamlines have been plotted for one quadrant in Figure 2. The detailed 
flow near the sphere is shown in Figure 3. Some interesting flow patterns have been 
predicted. While there is a thin layer on the surface of the particle, as in the case a 
particle at the velocity antinode, the streamlines are not closed. The reverse flows within 
the layer exist only for a belt-shaped region about the equator. Specifically, for angular 
positions in the range 55° < 9 < 125° with the polar axis, the thin layer of reverse flows 
is found. These values correspond to cos 2 9 = \ which represents the asymptotes of the 
stream function given by equation (37). Outside this range, around the polar regions, the 
streamlines merge with the exterior flow. Mathematically, this behavior results because 
the Stokes-layer solution consists of two components corresponding to the first and the 
third spherical harmonics. These solutions add up with two guaranteed zeros in the radial 
direction only for a limited region over the sphere as discussed. 

These results have some important implications for experimental studies relating to 
particle levitation. 

If we consider a sphere displaced from the antinode or the node of a standing wave, 
we need to expand the standing wave velocity u z such that 

u z = A cos[A:(2 + Zf))]e lT = [A cos kZo + A(k.z ) sin k.Zo + 0(k 2 z 2 ]e iT (55) 

for small kn, in which Z 0 is the displacement of the sphere from the antinode, and 
k = uj/c. Here, the first term of equation (55) is just the far-field velocity for Riley’s 
[9] problem, and the second term is the far-field velocity for the problem we have just 
discussed. That means the problem which the sphere is displaced from the antinode is 
the combination of the two problems we just discussed, together with some additional 
nonlinear terms. This is presently under investigation. 


References 

[1] R.E. Apfel. Technique for measuring the adiabatic compressibility, density, and 
sound speed of submicroliter liquid samples. J. Acoust. Soc. America, 59:339-343, 
1976. 

[2] A. Eller. Force on a bubble in a standing acoustic wave. J. Acoust. Soc. America. 
43:170-171, 1968. 

[3] A.R. Hanson, E.G. Domich, & H.S. Adams. Acoustic liquid droplet holder. Rev. 
Sci. Instrum., 35:1031-1034, 1964. 


14 



[4] L.V. King. On the acoustic radiation pressure on spheres. Proc. Roy. Soc. London , 
147:212-240, 1934. 

[5] L.V. King. On the theory of the inertia and diffraction corrections for the Rayleigh 
disc. Proc. Roy. Soc. London A, 153:17-40, 1935. 

[6] C.P. Lee & T.G. Wang. Outer acoustic streaming. J. Acous. Soc. America. 88:2367- 
2375, 1990. 

[7] RL. Marston & R.E. Apfel. Acoustically forced shape oscillations of hydrocarbon 
drops levitated in water. J. Colloid Interface Sci ., 68:280-286, 1979. 

[8] R.S. Rayleigh. Phil. Trans. Roy. Soc. London A, 245:535-, 1883. 

[9] N. Riley. On a sphere oscillating in a viscous fluid. Quart. J. Mech. Appl. Math., 
19:461-472, 1966. 

[10] H. Sclichting. Berechnung ebener periodischer Grenzschichtstrdmungen. Physik 
Zeitschr., XXXIII:327-335, 1932. 



Heat 

Transfer 

1994 

Proceedings of 

The Tenth International 
Heat Transfer Conference 

Brighton, UK 


Edited by 
G.F. Hewitt 

in cooperation with the members of the 
International Scientific Committee: 

M. Combamous, M. Cumo, E. Hahne, 

G. Hetsroni, C.J. Hoogendoom, 

J. R. Howell, A.I. Leontiev, J.S. Lee, 

L.F. Milanez, P. Oosthuizen, V.M.K. Sastri, 

K. Suzuki, B.X. Wang 


Volume 

5 

Pool boiling 

Particulates, porous media 
and special applications 
Natural and mixed 
convection 



HTD-Vol. 332 


Proceedings of the 

ASME Heat Transfer Division 

VOLUME 1 

• HEAT TRANSFER IN MICROGRAVITY SYSTEMS 

• RADIATIVE HEAT TRANSFER AND RADIATIVE HEAT TRANSFER 
IN LOW-TEMPERATURE ENVIRONMENTS 

• THERMAL CONTACT CONDUCTANCE AND 
INVERSE PROBLEMS IN HEAT TRANSFER 


presented at 

THE 1996 INTERNATIONAL MECHANICAL ENGINEERING CONGRESS AND EXPOSITION 
NOVEMBER 17-22, 1996 
ATLANTA, GEORGIA 


sponsored by 

THE HEAT TRANSFER DIVISION, ASME 


edited by 

ASHOK GOPINATH 

NAVAL POSTGRADUATE SCHOOL 

S. S. SADHAL 

UNIVERSITY OF SOUTHERN CALIFORNIA 

PETER D. JONES 
AUBURN UNIVERSITY 

J. SEYED-YAGOOBI 
TEXAS A&M UNIVERSITY 

KEITH A. WOODBURY 

THE UNIVERSITY OF ALABAMA 


THE AMERICAN SOCIETY OF MECHANICAL ENGINEERS 

United Engineering Center / 345 East 47th Street / New York, N.Y. 10017 



HEAT TRANSFER IN MICROGRAVITY SYSTEMS 


Introduction 
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Thermal phenomena in reduced gravity have attracted considerable scientific interest in a wide 
variety of topics including thermocapillary flows, crystal growth behavior, multiphase phenom- 
ena, containerless processing of drops and bubbles, and combustion, to name a few. This excit- 
ing branch of fluid mechanics and heat transfer offers a lot of promise for future scientific devel- 
opment, particularly in the creation of new types of material. The contributions from the heat 
transfer community have become an important element in this development. 

This symposium is part of a continuing series on Heat Transfer in Microgravity Systems spon- 
sored by the K-6 Committee (on Energy Systems) of the ASME Heat Transfer Division. The heat 
transfer community involved in microgravity research has been responding with considerable 
enthusiasm to our call for papers over the last few years. Ten papers were presented in two ses- 
sions — one primarily on phase-change phenomena, and the other mainly on thermocapillary 
flows and surface effects — at the 1996 ASME International Mechanical Engineering Congress 
and Exposition (Atlanta, Georgia, November 17-22, 1996). 

We are grateful to the authors who have participated in these sessions and offered their sci- 
entific contributions. We also thank the reviewers whose critical evaluations provided the neces- 
sary input for maintaining the technical quality of these symposia. 
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FOREWORD 


The availability of microgravity conditions has made possible certain thermal phenomena that 
were otherwise difficult to achieve. Among them is the deep undercooling of liquids, a state that 
is favorable for certain types of crystal growth. Such states may also be achieved by acoustic or 
electrostatic levitation, and studies involving containerless processing of materials are often as- 
sociated with microgravity studies. Low gravity also allows the isolation of various thermal phe- 
nomena that would otherwise be obscured by gravity driven flows. In particular, thermocapillary 
flows are a popular subject of investigation in this regard. Microgravity related research is currently 
active in a wide variety of analytical and experimental problems dealing with levitated drops and 
bubbles, thermocapillary flows, low gravity combustion, crystal growth, containerless processing, 
and thermophysical properties of undercooled liquids. 

This symposium volume consists of seven papers that were presented in two sessions. They 
cover several of the above topics, with participants from universities as well as government labo- 
ratories. We are thankful to all of the authors who have contributed to these sessions. We are also 
very grateful to the referees whose input has helped maintain the quality of ASME publications. 
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FOREWORD 


Thermal phenomena in reduced gravity have attracted considerable scientific interest, partic- 
ularly in deep undercooling and crystal growth. The availability of sustained microgravity condi- 
tions has made possible thermal phenomena that were otherwise difficult to achieve. In addition, 
low gravity experimentation has allowed the observation of certain thermal phenomena that 
have been obscured by the interference of gravity-dominated processes. Among them is the sub- 
ject of thermocapillary flows. Much of the work on undercooling has been possible with con- 
tainerless processing; this has also been achieved by acoustic or electrostatic levitation. Such 
ground based studies on containerless processing of materials are often placed in the realm of 
microgravity. Research activities related to microgravity presently cover various experimental 
and theoretical areas dealing with levitated drops and bubbles, containerless processing, ther- 
mophysical properties of undercooled liquids, thermocapillary flows, low gravity combustion, 
and crystal growth. 

This is the second in a series of symposia on microgravity sponsored by the K-6 Committee 
(Energy Systems) of the ASME Heat Transfer Division. Ten papers are presented in two sessions - one 
on general papers, the second on thermocapillary flows - at the 1994 International Mechanical 
Engineering Congress and Exposition (Chicago, Illinois, November 6-11, 1994). 

We are grateful to the authors who have participated in these sessions and offered their sci- 
entific contributions. We also thank the reviewers whose critical evaluations provided the neces- 
sary input for maintaining the technical quality of these symposia. 
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HEAT TRANSFER IN MICROGRAVITY SYSTEMS 


Introduction 


A. Gopinath 

Naval Postgraduate School 
Monterey, California 


S. S. Sadhal 

University of Southern California 
Los Angeles, California 


Low gravity research in heat transfer and fluid mechanics has always been of direct interest in 
space technology where the conditions necessitate the investigations of physical phenomena rel- 
evant to areas such as electronics cooling and combustion. More recently, there has been con- 
siderable scientific interest in relation to deep undercooling and crystal growth. With the avail- 
ability of sustained microgravity conditions, it has become possible to create certain thermal phe- 
nomena that were otherwise difficult to achieve. Among them is the formation of new types of 
crystals through deep undercooling in containerless systems. Also, experimentation under low 
gravity has allowed the observation of thermophysical phenomena, such as thermocapillary 
flows, that are normally obscured by gravity-dominated processes. Research activities related to 
microgravity are currently covering various experimental and theoretical areas dealing with lev- 
itated drops and bubbles, containerless processing, thermophysical properties of undercooled 
liquids, boiling phenomena, low gravity combustion, thermocapillary flows, and crystal growth. 

This symposium on microgravity is the third in a series sponsored by the K-6 Committee on 
Energy Systems of the ASME Heat Transfer Division. Eleven papers have been accepted for pub- 
lication, and these are being presented in two sessions. The first session consists of six papers 
on boiling phenomena in low gravity, including a review paper. In the second session there are 
three papers on other topics related to microgravity heat transfer. 

We appreciate the authors' participation through their scientific contributions to these ses- 
sions. At the same time, we are also grateful to the referees for providing critical reviews in order 
to maintain the technical quality of these contributions. 
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ABSTRACT 

Ultrasonic levitation allows the processing of low 
melting materials both in 1 G as well as in microgravity. 
The free suspension of the melts also facilitates 
undercooling, permitting the measurement of the physical 
properties of the metastable liquids. A convenient method 
to melt a levitated sample involves its spot heating 
through a focused radiant source, the heat input to the 
sample is controlled by the material emittance as well as 
the external convective flows. Because of high intensity 
sound fields required for levitation, thermoacoustic 
streaming will significantly increase the heat transfer from 
the sample to the environment, and it will therefore 
decrease the heating efficiency. Experimental measurement 
involving flow visualization and power input monitoring 
have allowed the quantitative assessment of this 
enhancement in heat transfer at ultrasonic frequencies and 
for millimeter-size samples. A decrease of temperature by 
up to 150 C for a sample initially at 550 C without the 
sound has been measured. Other results involving normal 
1 G and low gravity flow visualization and material 
processing are presented 


INTRODUCTION 

A low gravity environment is ideally suited for 
experimental studies involving the melting and 
solidification of materials in the absence of a container. 
The drastic reduction in the effects of the Earth 
gravitational acceleration allows the use of much lower 
levels in the electromagnetic (Okress et al., 1952, and 
Cummings and Blackburn, 1991), electrostatic (Rhim et 
al., 1993), or acoustic (Trinh, 1985) fields used to 
remotely position the sample of interest. Depending upon 


the material properties, the composition of the processing 
medium, and the levitation technique, the reduction of 
gravity could allow the study of phenomena not ordinarily 
observed, or it could lead to more accurate experimental 
measurements. For example, the decrease in the magnitude 
of the electromagnetic field required for positioning a 
molten metal droplet in low gravity could result in the 
enhanced supercooling of the liquid sample. The reduction 
of the acoustic intensity used for sample levitation in a 
gaseous host medium also means a reduced level of 
acoustically induced convective flow around the levitated 
sample. The effects of acoustic streaming associated with 
ultrasonic levitation is the subject of this experimental 
investigation. 

The specific results to be reported below have been 
obtained from a flow visualization study of streaming 
flow fields in a single-axis ultrasonic levitator operating at 
25 kHz and in 1 G. A quantitative assessment of the 
enhanced convective heat transfer from a locally heated 
sample has been obtained for moderately high sample 
temperature, and the results are summarized. Finally, 
some results of the levitation processing of low-melting 
materials in 1 G using ultrasonic techniques are described. 


I. THE EXPERIMENTAL TECHNIQUE 

A standard single-axis ultrasonic levitator (Trinh, 1985) 
is used to levitate the liquid and solid samples to be 
processed in a gaseous environment and at the focus of 
either a Nd-Yag laser or a Xenon arc lamp. Figure 1 
provides a schematic representation of the apparatus. An 
ultrasonic standing wave with an integer number of half- 
wavelengths (typically 1.5 or 3) is established between the 
driver and reflector, and a sample can be levitated at any 
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one of the pressure nodes (acoustic velocity antinodes). In 
addition to the second-order radiation pressure which is at 
the origin of the levitation force, a second order steady- 
state convective flow, or acoustic streaming field, is also 
generated. These flows are the ultrasonic equivalent of the 
streaming flow Fields investigated by Gopinath and Mills 
(1993) in the same context of the application of acoustic 
levitation. The visualization of these flows is carried out 
through light scattering by smoke particles under laser 
sheet illumination. The resulting streamlines are imaged 
by a video camera and recorded on tape. The temperature of 
the heated samples is determined through embedded 
thermocouples when they are mechanically suspended, and 
by the use of an infrared imaging camera when they are 
levitated. 



SCHEMATIC DESCRIPTION OF EXPERIMENTAL 
APPARATUS. THE SAMPLE (E) IS LEVITATED IN THE 
AXISYMMETRIC STANDING WAVE BETWEEN THE 
ULTRASONIC DRIVER (C) AND THE REFLECTOR <D). THE 
LIGHT SHEET (B) IS USED TO VISUALIZE THE FLOW 
FIELD THROUGH SCATTERING FROM SMOKE 
PARTICLES. A XENON ARC LAMP (A) BEAM IS 
FOCUSED ON THE SAMPLE THROUGH A LENS (F). 


II. FLOW VISUALIZATION OF 
THERMOACOUSTIC STREAMING 

Incense smoke was used as a tracer for these flow 
visualization studies. Previous investigations under 
isothermal conditions have revealed that the axial 
symmetry of the levitator was extended to the acoustic 
streaming flow field. This symmetry would not 
necessarily be preserved, however, in the case of a locally 


heated levitated or mechanically held sample because the 
spot heating is generally not along the axis of symmetry. 
With this caveat in mind, we have carried out flow 
visualization using sheet lighting parallel to the levitator 
axis of symmetry which is vertically oriented parallel to 
the gravity axis. The heated sample was simulated by a 
cylindrical mechanically held thermistor (1.9 mm diameter 
and 6.35 mm long). The temperature of the thermistor 
was varied by changing its drive voltage; its resistance 
was monitored by measuring the current input. All 
experiments were carried out at an ambient room 
temperature of about 23 C. The ultrasonic frequency was 
25 kHz and the maximum Sound Pressure Level or SPL 
was 155 dB. The SPL is defined as a logarithmic relative 
pressure measurement with a fixed reference pr e f = 0.0002 
liBar expressed in dB. 

Figure 2 is a photograph of a video frame showing the 
flow pattern around the unheated cylindrical thermistor 
with ultrasound at 145 dB. The sound driver and the 
reflector (separated by 1.5 acoustic half- wavelengths) can 
be seen at the bottom and top of the picture together with 
the primary and secondary eddies. An enclosure (6.35 x 
6.35 x 3.81 cm) has been placed around the levitation 
region in order to slow smoke dissipation. This enclosure 
does not significantly contribute to the primary standing 
wave in the region of interest around the sample. The 
primary set of eddies would be present in the chamber 
even without the thermistor, while the secondary eddies 
are attached to the thermistor and are caused by its 
presence. The primary eddy flow direction is counter- 
clockwise on the right side and clockwise on the left side 
of the picture; the direction of flow in the secondary eddies 
is directly opposite. Theory predicts a symmetrical set of 
vortices above and below the cylinder for streaming around 
a cylindrically shaped body. In practice, the superposition 
of the primary streaming flows caused by the enclosure 
reinforces the upper set of eddies, but opposes the lower 
vortices. 

Figure 3 shows a still video frame for a thermistor 
heated to 150 C in an ultrasonic standing wave with the 
SPL at 140 dB. In addition to the secondary vortices 
attached to the upper part of the thermistor, another set of 
vortices reappear on the lower side of the thermistor. This 
reappearance of the lower vortices is probably due to 
natural convective flows associated with the temperature 
gradient caused by the hot sample immersed in a cooler 
gas environment 

Figure 4 shows the streaming flow field immediately 
around the thermistor for a temperature of 450 C and the 
SPL at 145 dB. The upper eddies are no longer prominent, 
but the lower eddies are sharper and detached from the 
sample. These lower eddies are also split into two sets of 
counter-rotating vortices. The outer component of this 
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pair periodically stretches outward, sheds, and reforms 
closer to the inner component. The frequency of vortex 
shedding is controlled by the tuning of the ultrasonic 
standing wave. For a given sample temperature, a steady 
flow pattern can be established and maintained if careful 
tuning of the ultrasound is implemented. The influence of 
natural convection manifests itself through the shift in the 
location of the attached secondary vortices: free convection 
induces a flow opposite to the original primary streaming 
flow of the isothermal chamber and reinforces the lower 
set of secondary vortices while washing out the upper 
pair. 



FIGURE 2 

FLOW PATTERN AROUND THE UNHEATED CYLINDRICAL 
THERMISTOR WITH AN ULTRASONIC WAVE AT 145 dB. 




FIGURE 3 

FLOW PATTERN AROUND THE THERMISTOR HEATED FIGURE 4 

TO 150 C WITH SOUND AT 145 dB. FLOW PATTERN WITH THERMISTOR AT 450 C AND 

SOUND AT 145 dB 
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FIGURE 5 

TEMPERATURE OF A SPOT-HEATED CERAMIC SAMPLE WITH AND WITHOUT THE ULTRASOUND AT 160 dB 


III. THE EFFECT OF STREAMING ON FORCED 
CONVECTIVE HEAT TRANSFER 

Figure 5 displays the results of the measurement of a 
spot-heated ceramic sample mechanically held at the 
pressure node (velocity antinode) of an ultrasonic standing 
wave at 160 dB. The sample was heated by focusing the 
light from a Xenon arc lamp on one side of the sample. 
The temperature was obtained from a thermocouple 
embedded in the spherical sample of 3 mm diameter. A 
dramatic drop in the sample temperature from 640 to 440 
C (45% decrease) illustrates the effectiveness of forced 
cooling by ultrasonic streaming. Figure 6 shows that the 
effect of acoustic streaming is relatively large in the lower 
relative intensity range. Increasing the acoustic drive 
appears to indicate a "saturation" effect, as shown by the 
relatively smaller temperature drop at higher acoustic drive 
voltages. 

In order to simulate the microgravity conditions, 
additional measurements were carried out using the 
mechanically held thermistor and relatively low SPL (145 
dB). At 25 kHz and 145 dB in air at 25 C the Stokes 
boundary layer thickness 5 (S=(2v/io)°- 5 , where v is the 


kinematic viscosity of the fluid, and co is the angular 
frequency of the acoustic wave) is approximately equal to 
14 pm and the acoustic displacement (a = 2 jiU 1c /co) is 
about 35 Jim. U 1C is the acoustic particle velocity. 
Defining a streaming Reynolds number in the same 
manner as Gopinath and Mills (1993) and Leung et al. 
(1989): R, = (a/6) 2 , we find that the value for R, is 
approximately equal to 4 for our measurements at 143.5 
dB in air for a thermistor heated to 420 C. 

Assuming that the thermistor is isothermal, an average 
Nusselt number can be calculated. Nu = hd / k , where h 
is the average heat transfer coefficient, d is the thermistor 
diameter, and k is the thermal conductivity of air. The 
average heat transfer coefficient is based on a cylindrical 
geometry, and can be expressed as h = Q / \ (AT), where 
Q is the power input, A c is the surface area of the 
cylinder, and AT is the temperature difference between the 
thermistor and the ambient air. The results are listed in 
table I for four different temperature differences between 57 
and 397 C. Although the analysis of Gopinath and Mills 
is based on a spherical geometry, for small temperature 
differences, and for much larger values of R, , the ratio 
Nu / (R, J 0 - 3 was calculated for comparison. 
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Acoustic drive voltage (mV) 


FIGURE 6 

THERMISTOR TEMPERATURE AS A FUNCTION OF ACOUSTIC DRIVE FOR DIFFERENT INITIAL TEMPERATURES 


AT(°C) h(W/m 2 K) Nu=hd/K Nu/(R,) 0S 


57 

63 

4.13 

1.99 

132 

125.6 

7.02 

3.46 

307 

155.6 

6.70 

3.74 

397 

165.4 

6.41 

3.75 


TABLE I 

RESULTS OF THE CALCULATION OF THE AVERAGE 
HEAT TRANSFER COEFFICIENT AND THE AVERAGE 
NUSSELT NUMBER FOR FOUR DIFFERENT 
TEMPERATURES. 


IV. DISCUSSION 

Previous experimental studies of thermoacoustic 
streaming have also included flow visualization and heat 
transfer coefficient measurement (Fand and Kaye, 1961 , 
and Fand, 1962). Most of these results were obtained for a 
horizontally directed traveling sound wave normal to the 
axis of the cylinder and for a ratio of the acoustic half- 
wavelength to the cylinder diameter which is at least equal 
to 6 (X/2d is equal to 3 in this present work). The ratio of 
the acoustic particle displacement to the diameter of the 
cylinder a/d was equal to 0.016 in the previous work while 
it ranges from 0.019 to 0.028 in the current investigation. 
In addition the previous flow visualization results were 
obtained with flowing smoke streams while the present 
ones have been based on initially stationary smoke 
particles. The previous studies could not therefore reveal 
the superposition of isothermal streaming with the 
thermal component These differences explain why the 
current measurement of the average heat transfer 
coefficient few* acoustically forced convective cooling yield 
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results which are a factor of 5 larger than the values 
obtained in the earlier studies. 

The current flow visualization results suggest that 
natural convection plays a significant role in the 
determination of the morphology of the flow field: the 
resulting flow configuration is a superposition of the 
isothermal streaming field, natural convective flow, and 
the streaming flow associated with the sample itself. The 
principal controlling parameters are the temperature 
difference AT and the acoustic particle displacement 
amplitude relative to the sample diameter a/d. For 
moderate temperature differences ( AT < 500 C ), two 
distinct flow field regimes appear to dominate near the 
sample: a symmetric pattern of four eddies both above and 
below the sample at lower sound intensities (SPL < 145 
dB) and an asymmetric distribution where two sets of 
counter-rotating eddies are found at the lower half of the 
sample. Vortex shedding can happen, and it is probably 
related to the stability of the standing ultrasonic wave 
under thermal fluctuations. Careful tuning of the acoustic 
wave always appears to stabilize any large scale 
oscillation of the flow field. One must keep in mind that 
the quasi- isothermal streaming flow field characteristic of 
the chamber provides a background which must be taken 
into account in the detailed analysis of the overall fluid 
flow distribution. At very high sound levels ( SPL > 150 
dB) the smoke particle visualization technique is no longer 
effective because of the high velocities attained, and no 
clear structure can be resolved. It seems, however, that the 
lower vortex pattern attached to the sample still remains, 
even in the midst of what would appear to be turbulent 
flow. 


V. SAMPLE PROCESSING EXPERIMENTS 


Ultrasonic levitation melting, undercooling, and 
solidification in a gaseous environment of low melting 
metals, and organic or inorganic materials have been 
carried out in ground-based laboratories in isothermally 
heated chambers 9 . Spot heating of levitated samples is 
more difficult, however, because of the local thermal 
disturbance affecting an essentially rather well tuned 
system. The non-trivial task of maintaining levitation 
under the influence of acoustic streaming has been 
investigated in 1 G, and it can be accomplished for a 
sample-environment temperature difference of about 500 C 
for low density materials. This implies that the 
microgravity implementation of such a capability has 
been verified. The extension to higher temperature 
processing could be accomplished through the 
combination of a high temperature isothermal facility 
coupled to a spot heating capability. 


The spot heating of materials for containerless 
processing is advantageous in terms of high heating and 
cooling rates, but it presents the disadvantage of requiring 
the implementation of reliable remote temperature 
measurement capabilities. In addition to standard single 
and multi-color pyrometers, we have chosen to implement 
infra-red imaging cameras for high temperature Earth-based 
developmental studies as well as for experiments in 
material processing. 

Figure 7 reproduces video displays of an infra-red 
imaging camera recorded during the levitation melting of a 
spherical polymer sample in 1 G. The sample starts to 
soften and increases in volume due to boiling of volatile 
components. It ultimately distorts and solidifies into a 
non spherical sample. 

Because of the thermal input from the spot heating, the 
ultrasonic standing wave undergoes oscillations which 
cause sample translational vibrations. These positional 
instabilities (also accompanied by rotational motion) in 
turn initiate thermal fluctuations because the sample has a 
time-dependent motion into and out of the heating beam 
or focal region. Earth-based processing is more sensitive 
to this instability process because of the high sound 
pressure levels required for levitation. The practical 
solution to this problem involves real-time, fast response 
ultrasound retuning as well as slow increases in the heat 
input A closed-loop feedback system coupling the sound 
tuning and intensity to the sample motion and heat input 
is required for Earth-based controlled processing. The 
sample behavior in microgravity should remain 
qualitatively similar except for slower position 
fluctuations due to the smaller magnitude restoring force. 
The same active control capabilities for the levitator and 
heat input, however, will still be required. 

Figure 8 are video still frames of the IR imager output 
for the Earth-based levitation heating of a low density 
shuttle tile sample (3 mm in diameter). The series of 
photographs shows the typical fluctuations in position of 
the sample. The maximum temperature recorded is 
approximately 520 to 550 C. A substantial thermal 
gradient exists on the surface of the sample as it is heated 
from one direction by a 2 mm diameter gaussian beam 
from a Nd-YAG CW laser. 

Figure 9 shows a video still frame sequence describing 
the various stages of melting and solidification of an O- 
Terphenyl sample in low gravity during the parabolic 
flight of the NASA KC-135 airplane. About 15 seconds 
are available at an acceleration level of about 0.05 G for 
positioning, melting, and solidifying the sample. The 
resulting oblate shape is caused by the higher intensity 
sound required for levitation under 1.8 G during the high 
G part of the parabolic flighL 
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These preliminary experiments have provided strong 
evidence for the significant role played by acoustic 
streaming in the levitation processing of materials. Flow 
visualization has shown that the resulting convective field 
is a combination of forced isothermal streaming 
characteristic of the chamber, near sample boundary 
streaming, and finally natural buoyancy flows. The 
reduction of the gravitational component thus 
significantly alters the flow morphology by virtually 
eliminating the natural buoyancy component and by 
reducing the maximum acoustic field intensity required for 
sample positioning. 

Results of the measurements carried out in this work 
also suggest that the enhancement of heat transfer quickly 
increases for sound pressure levels up to 150 dB, but 
levels out at the very high SPL. Since the lower end of 
the sound intensity scale is likely to be implemented in 
low gravity, the effect of thermo-acoustic streaming 
should remain an important factor in microgravity 
processing of materials in gases. Acoustic streaming 
affects not only the efficiency of spot heating of freely 
suspended samples, but it also significandy degrades their 
position and rotational stability within the ultrasonic 
field. An appropriately designed feedback system tying the 
sound tuning and intensity and heat input controls to 
sample position information will be required for both 
ground-based and microgravity processing. 

Earth-based processing of low-melting polymeric and 
metallic materials has provided early evidence of the 
feasibility of spot heating for moderate sample- 
environment temperature difference. An experimentally 
determined upper limit of about 500 C for the sample- 
environment temperature difference has been obtained for 
ground-based processing of low temperature samples. The 
environment temperature can of course by as high as 
required as long as isothermal conditions are maintained 
and the positioning (levitation) system remains effective. 


FIGURE 7 


LEVITATION PROCESSING OF A SPHERICAL 
POLYMERIC SAMPLE IN AN ULTRASONIC FACILITY 
WITH FOCUSED XENON ARC LAMP HEATING 
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FIGURE 8 

IR IMAGER OUTPUT DURING THE SPOT HEATING OF A LEVITATED SHUTTLE TILE SAMPLE IN 1 G. 
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LOW GRAVITY MELTING AND SOLIDIFICATION OF O-TERPHENYL DURING PARABOLIC AIRPLANE FLIGHT 
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ABSTRACT 

The thermal analysis of the effect of the steady 
streaming motion induced by a solid sphere in an acous- 
tic field is carried out by the singular perturbation 
method. For sufficiently large frequencies, a thin Stokes 
layer region on the surface of the sphere constitutes the 
inner solution for the flow field which is then matched 
with a suitable outer solution. The temperature field 
in the fluid is induced by a thermal oscillation on the 
sphere boundary while the ambient is held at a con- 
stant temperature. The interaction of the thermal os- 
cillations and the acoustic field results in a nonzero 
time averaged steady convective transport of heat in 
the fluid. This is the principal effect being investi- 
gated in this study, with attention restricted to gases 
with Pr % 1. As in the case of the velocity field, a 
singular perturbation procedure is employed for deter- 
mining the temperature distribution in the fluid. For 
large streaming Reynolds numbers, R s , being consid- 
ered, the matching of the inner Stokes layer with the 
outer field is through a thicker outer boundary layer. 
While the inner solution can be found analytically, the 
outer solution requires numerical work. The final re- 
sults show that although there is no net exchange of 
heat between the sphere and the fluid, there is, how- 
ever, a steady flow of heat, into and out of, each hemi- 
sphere of the physical domain. This flow causes an 
equal heating and cooling of the fluid in each hemi- 
sphere, while within the solid sphere it is indicative of 
a steady temperature gradient across its poles. 

NOMENCLATURE 

a Radius of the sphere 

A Oscillation amplitude of the sound wave, Uoo/u 

c Velocity of propagation of sound in the fluid 
A/ 2 Frequency parameter, a 2 jj/u 

N iii Inner average Nusselt number 

Nuq Outer average Nusselt number 

p m Mean fluid pressure 


po Acoustic pressure amplitude 
Pr Prandtl number 
r Radial coordinate, r m /a 
R s Streaming Reynolds number, = s 2 M 2 

t Dimensionless temperature defined in §4 
T Dimensional temperature of the fluid 
Temperature of the ambient 
(AT) a Amplitude of temperature oscillations 
u,v “Artificial” velocities defined in §4 
u r Angular, radial velocity components 
Uoo Velocity amplitude in the sound wave, Alj 
x, y Boundary layer coordinates defined in §4 
z * Coordinate in the axial direction 

Greek Letters 

7 Phase difference between oscillations 
7o, 71 Phase angles defined in Eqs. 19 and 20 
7r Ratio of specific heats of the fluid 
6 Stokes layer thickness, y f vju 
e Amplitude parameter, A/a 
r], i) Normal boundary layer coordinates 
p, 9 Angular coordinates, p = cos# 

A Wavelength of the sound field 
v Kinematic viscosity of the fluid 
r Dimensionless time, wr' 

<t> y $ Dimensionless temp, diff., (T — Too)/(AT) a 
0, V Dimensionless stream function 
w Angular frequency 

Superscripts 

* Represents dimensional quantities 


1. INTRODUCTION 

In this study, we examine the process of convective 
heat transfer due to acoustic streaming induced by a 
sound field about an isolated sphere which is subject 
to time-periodic temperature fluctuations. The prin- 
cipal feature of interest in the present study is the en- 
ergy transport phenomenon emanating from the time- 
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independent contribution of the convective term due 
to the interactions of the thermal oscillations with the 
acoustic field. 

In previous studies (Gopinath & Mills, 1993, here- 
after referred to as [I] , and Gopinath, 1993) the steady 
heat transport due to the streaming motion was exam- 
ined for the case of an isothermal sphere exchanging 
heat with an isothermal fluid. However, in the cur- 
rent study, although the fluid ambient is considered 
isothermal, the temperature of the sphere is taken to 
be time-periodic. Such a situation could arise in, say, a 
pulsed heating process (see Magill et al., 1987) or other 
inherently time-periodic heat transfer processes. The 
basic problem of the sphere treated here would help ad- 
dress the larger issue of the influence of acoustic fields 
on such processes. Unlike the work of Arpaci (1966, 
pp. 324-335), or more recently Khedari et al. (1992), 
who have considered the conductive transfer of heat 
in solid bodies subject to time-periodic temperature 
boundary conditions, the present work deals princi- 
pally with the convective transport of heat due to steady 
thermoacoustic streaming effects in the fluid around 
the sphere. 

The steady streaming motion in the fluid is taken 
to be induced around a rigid sphere of radius, a, by a 
standing acoustic field with a velocity distribution of 
the fortn, Uoo cos(u»r*) sin(2?rz*/A). The development 
follows that of [I] wherein length scale arguments were 
used to treat the parameter ranges of au/c <3C 1 and 
e = (Zoo /aw «: 1, for which the flow induced around 
the sphere can be assumed to be laminar and unsep- 
arated, with negligible compressibility effects. This 
also allows the governing equations to be treated by 
the method of matched asymptotic expansions with 
e playing the role of a small perturbation parameter. 
Furthermore, only the high frequency range, M 2 = 
a 2 (jj/u ^ 1, is considered, for which the streaming ef- 
fects are most significant. For the above conditions, 
the basic flow field around the sphere was developed 
by Riley (1966) who showed that the steady flow field 
of interest can be divided into essentially two distinct 
regions — a thin, inner recirculating Stokes layer re- 
gion in which an O(eUoo) streaming velocity originates 
to drive the steady flow in an outer region making up 
of the remainder of the domain. Stuart (1966) recog- 
nized that the strength of the steady flow in the outer 
region is governed by a streaming Reynolds number, 
R $ = £ 2 M 2 , which can be uniquely determined from 
the driving acoustic signal (see [I]). In the present 
study, the acoustic signal is taken to be sufficiently 
strong so as to give rise to large R $ , for which the 
steady transport effects due to the streaming motion 
are most pronounced. Particular attention -is given 
to cases in which the surrounding fluid is a gas with 


Pr = 0(1). For large R 9y the outer steady flow has a 
boundary layer structure (Stuart, 1966), the behavior 
of which has been obtained in [I] . 

For purposes of simplicity the periodic tempera- 
ture excursions of the sphere are taken to be har- 
monic (at a single frequency) and of the form, + 
(A T) a cos(wr* +7). It is assumed that the amplitude 
of these oscillations, (AT) a , is small enough to ne- 
glect (as a first approximation) any interaction of the 
thermal and acoustic fields. Furthermore, it is also as- 
sumed that any high-intensity thermoacoustic effects 
as discussed by Gopinath (1993) are small. The va- 
lidity of this assumption can be ensured if a suitably 
defined Eckert number is maintained small. More im- 
portantly, it must be noted that the angular frequency 
of the acoustic field, w, is taken to be “tuned” to match 
that of the temperature oscillations, with allowance 
made for a possible difference in phase, 7. Such iso- 
harmonic situations are of principal interest since the 
magnitude of the expected steady heat transport re- 
sulting from the interactions of these oscillations is the 
strongest for such cases. In general however, the pro- 
cedure followed in this study can be extended to any 
arbitrary periodic temperature disturbance, after it is 
Fourier decomposed into its constituent frequencies. 


2. GOVERNING EQUATIONS 


For this axisymmetric problem, we describe the 
fluid motion by the Stokes stream function (\p = ip* / 
Uood 2 ) in a spherical coordinate system, i.e., 


1 dti> 

Ur = 2^“ and 

r 2 o/j 


u* = - 


(l-/j 2 ) * di> 


dr' 


( 1 ) 


where /i = cos0. The dimensionless governing equa- 
tions of momentum and energy for the fluid are, 

Tr^*y + P K + ! 2£)1 H = W 13 '*' 

( 2 ) 


and 


37 + ^ [i$$\ ■ prm* 7 *' (3) 

where the operators D 2 , L and V 2 are defined as 


D 2 = 


d 2 (1 - / j 2 ) d 2 


dr 2 


dfi 2 ’ 


L = 


d 1 d 
■ + -■ 


(1 — /i 2 ) dr r dn ' 

(4) 


and 

2 1 d r 2 d 1 Id 

V ~ r 2 dr [ r dr] + r 2 dp 




djt 


(5) 


The boundary conditions are 


dip 


<f> = cos(r + 7 ) and ip — = 0 at r = 1 ( 6 ) 
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0 -4 ±r 2 (l - /j 2 )cost 

<j> -► 0 

In this study we are concerned with the cases of 5 < 1 
and M > 1. This parameter range has been treated 
by Riley (1966) to develop a solution for the flow- 
field in terms of inner and outer complementary per- 
turbation series expansions in powers of e and 1/M. 
It has been further supplemented in [I] with a treat- 
ment of the large R s flow field and the associated heat 
transfer problem. Similar methods are used in the fol- 
lowing sections with a focus on determining the ad- 
ditional steady thermoacoustic streaming phenomena 
introduced by the interaction of the oscillating thermal 
and acoustic fields. These methods call for perturba- 
tion expansions in the inner and outer regions. 

3. THE INNER REGION 


as r 


‘00 


(7) 


for which the leading order contribution was obtained 
by Riley (1966, Eq. 28) as, 


3\/2 o\ 

*0 = 


x \y2 77COST - cos(r — + e n cos(r - rj - (13) 


Although the above form for ♦o will suffice for the 
present study, it is important to note that it is the 
time-independent part of the 0(f) contribution, 
(Riley, 1966, Eq. 32) which explains the behavior of 
the acoustic streaming motion in the inner region and 
provides a description of the slip— like velocity which 
drives the steady flow in the outer region. 

For the temperature distribution also, a perturba- 
tion expansion of the form 


<£ = $o + 5($i, + *i„) + ■ ■ • (14) 


For large values of the frequency parameter, M, 
the basic oscillatory flow has a Stokes boundary layer 
on the surface of the sphere with a largely irrotational 
exterior region. The mechanics of this inner oscillatory 
shear layer region of dimensional thickness of 0{6) is 
best understood in terms of suitably defined Stokes 
boundary layer variables (Riley, 1966), 


\f 


\f 


r7=(r-l)-^= and ^(rj, /i, r) = V f (r, /i, r) (8) 


with <f>(7j,/i,r) = <2>(r,y/,r), in view of which the gov- 
erning Eqs. (2)-(3) become, 


d_ 

dr 


/<9 2 tf\ f d(*<d 2 */dri 2 ) 2 /j <9* d 2 *' 

V"(V/ L d(r},n) + 1 - n 2 dr) dvr . 

= \ $T + OUA/-‘,Ar 3 ) (9) 

2 or/* 


dr 


+ s 


>) ' 

. d(ri,n) _ 


1 


d 2 $ 


+o(fA/ _i : 

2\/2 d<P 


2 • Pr 


dr? M 


+ 0{ A/- 2 ) 


( 10 ) 


The above equations can be made to satisfy only the 
inner boundary conditions on the sphere surface, namely, 


is sought. Here it can be shown from a solution of 
Eq. (3) subject to Eqs. (6)— (7) that the leading order 
solution, 4> 0 , is given by 

<l>o = e~' ,v ^cos(r + 7 - r/\/Pr) (15) 

and represents an oscillatory temperature wave in the 
Stokes layer region. Since the variation of $0 is time- 
periodic, there is no net time averaged transfer of heat 
at this level. Of greater interest is the O(e) contribu- 
tion to the temperature distribution in the fluid which 
arises only due to the presence of the acoustic field. An 
analysis of the 0(e) balance of Eq. (10) using Eq. (14) 
shows that in addition to an 0(1) second harmonic and 
an 0{\ls/Rl) first harmonic contained in <&i u , there 
also arises from the convective term a non zero steady 
part, <ti,, resulting from the time-averaged interaction 
of the leading order oscillatory temperature and flow 
fields in the fluid. The behavior of this steady part is 
governed by 

1 / fl(*o,*o) \ (16) 

2 Pr dr) 2 \ d{i),n) r 

where the angle-brackets, ( ), denote a time-average 
of the enclosed quantities. This equation is subject to 
the inner boundary condition 

*„=0 at r; = 0 (17) 


<l» = — = 0 and <!> = cos(r + 7) at 17 = 0 
an 

( 11 ) 

The solution for the stream function, ♦, is sought in a 
perturbation expansion of the form, 

= ♦<, + + *iu) + •• • (72) 


and appropriate matching with the outer region. Al- 
though the details of the analysis have been excluded, 
it can be shown that the resulting variation of $1, can 
be obtained from a solution of the above equations as, 

(J> u = ^ g-i'/P? ^yp7(l - r7)cos(7 - 77\/Pr) 

- (2 + ij\/Pr) sin(y - q\/ Pr) 
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+ , 1 PrV ^ 2 e " {(1 -Pr)cos(7 + 7;-r ? v / Pr) 

(1 + Pr) v 

+ 2v , 'Pr sin(7 + 7 — ij\/Pr))J 
+ M*ioo(Pr,7) (18) 


where 


*ioo(Pr,7) 


with tan 70 


3 x/9Pr + 4 


2 (1 + Pr) 
n/Pt(1 + 3Pr) 
2(1 + 2Pr) 


sin(7 - 70) 


(19) 


There are two important implications of this steady 
temperature variation, both of which depend strongly 
on the phase difference, 7, namely : 


(1) Equation (18) predicts a non-zero fluid temper- 
ature gradient at the sphere wall given by 


d*u 

dv 


= sin(7i + jt/ 4) sin(7 - 71) 

rjzzO V 2 


with 


_ n/Pt - 1 

+ 1 


(20) 


This yields an average inner Nusselt number (bas- 
ed on the sphere diameter) for the corresponding 
heat transfer rate over each hemisphere as, 


Nui 

\ fR, 



-- sin(7i + >r/4) sin(7 - 71) (21) 


(2) The variation in Eq. (19) also prescribes a tem- 
perature at the outer edge of the inner Stokes 
region (*u(^-+oo)->^i $\oo) which in turn de- 
termines the steady temperature distribution in 
the outer region. 


4. THE OUTER REGION, R, » 1 

For the related problem of the steady variation of 
velocity and temperature in the outer region, the so- 
lution in the form of series expansions similar to Eqs. 
(12) and (14) is sought. Therefore, the following ex- 
pansions are used 

0 = 0o + £ (0u + 0i«) + . . • (22) 

0 = 0o + £ (0i# + 0iu) + ■ • • (23) 

with particular interest in the dominant contributions 
of the time-independent portions of each, namely the 


steady component of the 0(e) term, 0i # , in Eq. (22) 
for the stream function and the 0(e) term, 0i,, in 
Eq. (23) for the temperature. For completeness it is 
also useful to quote the solution for the basic lead- 
ing order contribution to the stream function, 0o in 
Eq. (22), obtained by Riley (1966, Eq. 27) as, 

^1 - -^) (1 -/j 2 )cosr (24) 

For the leading order temperature, 0o in Eq. (23), 
matching (as r/— ► oo) with the exponentially decaying 
oscillatory behavior of the inner region in Eq. (15) 
shows that in the outer region 0o = 0. 

Stuart (1966) showed that the nature of the steady 
transport effects in the outer region is governed by 
the magnitude of the streaming Reynolds number, R t . 
This parameter can be determined from the acoustic 
signal and the fluid properties as described (for the 
case of air) in Eqs. ( 13)— ( 16) of [I] . In general, for a 
plane standing acoustic field in an ideal gas, R, can be 
expressed as 



where (po/Pm) is the pressure amplitude ratio, which 
is oftentimes the measured parameter used to charac- 
terize strong acoustic fields. 

For large values of R s and Pr = 0(1) being con- 
sidered here, the variation of the steady terms \p u and 
01, in this outer region exhibits a boundary layer be- 
havior (Stuart 1966), as mentioned before in §1. Al- 
though this boundary layer region is thin on the scale 
of the sphere radius, it is much thicker than the in- 
ner Stokes layer. This region has to be analyzed by 
a numertcal solution of the governing equations sub- 
ject to suitable matching conditions from the inner re- 
gion. These equations themselves may be obtained in 
a manner described by Riley (1966). The process is 
not entirely trivial and the details have been omitted 
here. The relevant equations as developed in [I] have 
been summarized below in their final forms. The outer 
boundary layer variables are defined as, 

i) = (r-l)s/W, (26) 


and 

tf'i. (»>./*) = 

(27) 

For convenience these are expressed in terms of com- 
monly used symbols for the coordinates and suitably 
defined “artifical" velocities, 


y = V . x = fi,u = 



di>i, 


(28) 



along with 

t(i, y) = <j>i,(n, fj). (29) 

The governing equations in these variables are, 


du dv 
dx + dy 


du du 

Tt + "Ty + 


0 

(T^T 


dt dt 
U Yx +V Ty 


with boundary conditions 


0 (30) 


&u 

dy 2 

i a 2 < 

Prdy 2 


(31) 

(32) 


u = x (1 — x 2 ), v = 0, t — j $ loo at y — 0 (33) 

u — ► 0, t -> 0 as y oc (34) 

It may be recalled from the brief discussion at the end 
of §3, that the limiting value of the temperature from 
the inner region, $ioo in Eq. (19), prescribes the driv- 
ing temperature for the outer region as is clear from 
the temperature boundary condition in Eq. (33). 

The governing equations are now completely de- 
fined and can be solved with the help of a suitable nu- 
merical method. The coupled set of nonlinear partial 
differential equations for u and v are first solved using 
an implicit finite-difference scheme with marching of 
the solution from the equator (x = 0) to the poles (x = 
±1). All the derivatives are approximated by central 
differences and the nonlinearity is handled by quasi- 
linearization and iteration at each x-station along the 
periphery. Owing to the decoupled nature of the mo- 
mentum and energy equations, once u and v have been 
determined, the temperature, can be found in a rela- 
tively straightforward non-iterative manner using tri- 
angular resolution and backward substitution. The 
svmmetry about the equatorial plane is exploited to 
carry out the procedure over only one hemisphere. 
Further details of the discretization scheme and the 
grid parameters as well as an account of the means of 
accomodating the converging flow in the polar regions 
may be found in [I] . 

It must be emphasized that the system of gov- 
erning equations and boundary conditions presented 
above for the outer region, is strictly valid only for 
cases of strong streaming motion (R s » 1) in moder- 
ate Prandtl number fluids (Pr ~~ 1). 


5. RESULTS AND DISCUSSION 

It should be noted at the outset that the results 
for the steady flow field, are known and a repre- 
sentative plot of the boundary layer velocity profiles 
may be found in Fig. 3 of [I] . It is the variation of 


the steady temperature and the related heat transfer 
effects which will be of special interest in this study. 

A representative plot of the numerically determined 
outer boundary layer temperature profiles for air (Pr = 
0.7) is given in Fig. 1 for the case of (7 — 70) = ^/2 and 
0° < 9 < 90°. For the other hemisphere (90° < 9 < 
180°), the sign on the temperature values is reversed 
and the corresponding distribution may be simply ob- 
tained by a reflection about the t/-axis. 

Also of interest, is the variation of the local heat 
flux over the periphery of the sphere. This is character- 
ized by the local driving temperature gradient for the 
outer region, (— dt/dy) y =o, which is plotted in Fig. 2, 
also for Pr = 0.7 and (7 - 7o) = */2. The observed 
trend in this figure may be explained as follows: since, 
according to Eq. (33), / - r, the magnitude of the 
temperature and its gradient increase with |x| as the 
flow progresses from the equator to the poles. How- 
ever, as the flow converges towards the poles (x = ±1), 
continuity dictates a thickening of the boundary lay- 
ers. This in turn reduces the driving fluid temperature 
gradient, which becomes negligible as 1 -> ±1, where 
there is a breakdown of the boundary layer structure 
of the flow. The temperature gradient thus reaches 
a maximum at some intermediate angular location, 
which for air (Pr = 0.7), occurs at 9 « 45°, 135°. The 
resulting heat transfer rate with the fluid (over each 
hemisphere) can be characterized by an average outer 
Nusselt number (based on the sphere diameter) by nu- 
merically integrating this driving temperature gradi- 
ent. For air (Pr = 0.7) this gives 


jV — = f dx « 1.20 sin(7 - 70) (35) 

z\fW s Jo V ^y/y= 0 

An observation of the /i-dependence of the tem- 
perature in Eqs. (18)— (20) (and hence in Eq. (33)) 
shows that the driving temperature and its gradient 
in both the inner and outer regions, are antisymmetric 
about the plane of the equator. This indicates that 
the Nusselt number results in Eqs. (21) and (35) are 
only valid for each hemisphere, and there is no net 
exchange of heat between the entire sphere and the 
fluid. Such a situation is physically realized in the 
fluid by an equal amount of cooling and heating in each 
hemispherical portion of the domain, while within the 
solid sphere it takes the form of a steady flow of heat 
into a hemisphere and out the other, across the equa- 
torial plane. Thus in the sphere, this time averaged 
heat flow rate is capable of inducing a steady temper- 
ature gradient across its poles. It is also clear from 
these Nusselt number results that the magnitude and 
direction of this heat flow strongly depend on the re- 
lation between the imposed phase difference, 7, and 
the innate phasing provided by the fluid via 70 and 
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71, which is a characteristic feature of such thermoa- 
coustic flows. We are currently investigating the roles 
of the solid/lluid thermal properties in such a ther- 
moacoustic interaction, for a physically more realis- 
tic situation where a controlled pulsating surface heat 
source (rather than a temperature source) as discussed 
by Magill et al. (1987) is applied. 

In concluding, it is emphasized that this fundamen- 
tal problem serves to underline the importance of the 
ability to induce steady heat transfer rates to/from 
a body subject to time-periodic temperature fluctua- 
tions. A suitably chosen acoustic field is an important 
participant here and without it, the body does not ex- 
perience a steady exchange of heat with the surround- 
ing fluid. 
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Figure 1. Outer boundary layer temperature profile* 
for air (Pr = 0.7) for the case of (7 - 70 ) = */2 and 
large R (The corresponding angular locations, 0, in 
degrees are shown on the figure) 



Figure 2* Distribution of the local outer heat flux for 
air (Pr = 0.7) over the upper hemisphere for the case 
of (7 - 70) = ?r/2 and large R». 
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